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Preface 


This  report  describes  results  and  some  consequences  of  classifying  shallow- 
water  wind  wave  energy  spectra  in  terms  of  characteristic  wave  height,  fre¬ 
quency,  and  direction.  The  work  was  motivated  by  the  need  to  quantify  the 
variability  in  natural  energy  distributions  when  only  three  parameters  are  used 
to  define  wave  conditions.  A  special,  high-resolution  directional  wave  gauge 
is  used  for  this  purpose.  This  effort  was  authorized  by  Headquarters,  U.S. 
Army  Corps  of  Engineers  (HQUSACE),  under  Civil  Works  Coastal  Flooding 
and  Storm  Protection  Program  Research  Work  Unit  32484,  "Directionality  of 
Waves  in  Shallow  Water."  Funds  were  provided  through  the  Coastal  Engi¬ 
neering  Research  Center  (CERC),  U.S.  Army  Engineer  Waterways  Experi¬ 
ment  Station  (WES),  under  the  program  management  of  Ms.  Carolyn  M. 
Holmes,  CERC.  Messrs.  John  H.  Lockhart,  Jr.,  John  G.  Housley,  Barry  W. 
Holliday,  and  David  A.  Roellig  were  HQUSACE  Technical  Monitors. 

This  report  was  prepared  by  Dr.  Charles  E.  Long  at  CERC’s  Field  Re¬ 
search  Facility  (FRF)  in  Duck,  NC,  under  the  direct  supervision  of  Mr.  Wil¬ 
liam  A.  Birkemeier,  Chief,  FRF,  and  Mr.  Thomas  W.  Richardson,  Chief, 
Engineering  Development  Division,  CERC.  General  supervision  was  provid¬ 
ed  by  Dr.  James  R.  Houston  and  Mr.  Charles  C.  Calhoun,  Jr.,  Director  and 
Assistant  Director,  CERC,  respectively. 

The  directional  wave  gauge  and  its  data  processing  software  were  designed 
by  Dr.  Joan  M.  Oltman-Shay  while  at  Oregon  State  University  working 
through  an  Intergovemmentd  Personnel  Agreement.  Physical  maintenance  of 
the  gauge  was  provided  by  the  FRF  dive  team  consisting  of  Messrs. 
Birkemeier,  Michael  W.  Leffler,  H.  Carl  Miller,  Eugene  W.  Bichner,  and 
Brian  L.  Scarborough.  Gauge  calibration  was  maintained  by  Mr.  Kent  K. 
Hathaway,  FRF.  Data  acquisition  and  storage  were  performed  by  Mr. 

Clifford  F.  Baron,  FRF.  This  document  was  edited  by  Ms.  Janean  Shirley, 
Information  Technology  Laboratory,  WES.  The  contributions  of  all  of  these 
individuals  are  gratefully  acknowl^ged. 

At  the  time  of  publication  of  this  report,  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


1  Introduction 


Descriptions  of  ocean  wind  wave  conditions  in  coastal  engineering  studies 
are  often  restricted,  either  by  common  practice  or  lack  of  instrumentation,  to 
time  series  of  three  characteristic  parameters:  a  wave  height,  a  frequency 
associated  with  the  most  energetic  waves,  and  a  propagation  direction  that  is 
typical  of  the  most  energetic  waves.  Knowledge  of  these  three  parameters  for 
a  given  time  is  sufficient  to  employ  a  great  deal  of  the  guidance  in  the  Shore 
Protection  Manual  (SPM)  (1984),  perhaps  the  most  widely  referenced  treatise 
on  coastal  engineering  practice.  In  cases  where  wind  wave  energy  is  concen¬ 
trated  in  a  narrow  band  of  frequencies  and  a  narrow  range  of  directions,  a 
three-parameter  characterization  of  sea  state  is  reasonable,  and  engineering 
computations  based  on  SPM  guidance  may  be  considered  fair  estimates  of 
natural  processes.  However,  it  is  known  from  at  least  one  detailed  study  of 
shallow-water  wind  waves  that  energy  is  commonly  spread  over  broad  ranges 
of  frequency  and  direction  (Long  and  Oltman-Shay  1991).  In  these  cases,  a 
three-parameter  definition  of  wave  conditions  is  likely  to  lead  to  considerable 
errors  in  engineering  estimates.  The  intent  of  this  report  is  to  elucidate  and 
quantify  some  of  the  possible  errors. 

Quantification  is  achieved  from  high-resolution  empirical  observations  of 
wind  wave  frequency-direction  energy  spectra  obtained  in  shallow  water.  A 
set  of  6,759  spectra  from  regular  observations  over  a  5-year  time  span  is 
classified  and  grouped  by  ranges  of  characteristic  wave  height,  frequency,  and 
direction.  The  number  of  constituent  spectra  in  these  groups  illustrates  the 
variety  of  spectra  in  the  database  and  provides  some  measure  of  their  frequen¬ 
cy  of  occurrence.  Properties  of  the  constituent  spectra  in  each  group  are 
compared  to  corresponding  properties  of  a  unimodal,  monochromatic  wave 
train,  which  is  the  fundamental  wave  field  defined  by  just  one  height,  frequen¬ 
cy,  and  direction.  Spectral  properties  include  mean  spectral  shape,  a  measure 
of  variability  of  spectral  shape,  estimates  of  longshore  energy  flux,  and  esti¬ 
mates  of  radiation  stress  tensor  components.  All  of  these  properties  are  ger¬ 
mane  to  a  definition  of  nearshore  wave  climate  and  to  fundamental  wave 
properties  used  to  model  wave  forces  on  structures,  wave-generated  currents, 
and  wave-induced  sediment  transport. 

Because  this  study  is  based  on  empirical  observations  from  a  single  site,  it 
is  emphasized  strongly  at  the  outset  that  all  quantification  is  specific  to  that 
single  site.  Conclusions  reached  here  can  be  applied  qualitatively  to  other 
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sites  having  comparable  geometries  and  boundary  conditions,  but  should  not 
be  oonstni^  to  be  universally  ^plicable  to  all  sites. 

This  report  is  organized  into  six  chapters.  Chapter  2  describes  the  experi- 
mem  site,  the  directional  wave  gauge,  the  spectral  classification  scheme,  and 
the  distribution  of  constituent  spectra  among  the  classes.  Chapter  3  describes 
characteristic  spectral  shtq[>es  and  their  variability  within  classes.  Chapter  4 
shows  comparisons  of  estimates  of  longshore  energy  flux  from  constituent 
directional  spectra  and  from  corresponding  unidirectional  wave  trains  for  each 
of  the  classes.  Chapter  5  does  a  similar  analysis  for  estimates  of  components 
of  the  radiation  stress  tensor.  Chapters  4  and  S  both  discuss  quantifiable  bias 
errors  and  the  magnitudes  of  residual  errors  that  occur  when  a  three-parameter 
definition  of  sea  state  is  used.  Ch2q>ter  6  summarizes  salient  results  and 
makes  suggestions  for  future  work. 

Though  this  report  emphasizes  differences  in  certain  computations  when 
two  forms  of  sea  state  description  (complete  spectrum  versus  three  parame¬ 
ters)  are  used,  it  also  provides  one  form  of  reduaion  of  a  rather  large  spectral 
database  into  meaninghil  classes  of  spectral  shapes.  These  shapes  represent 
eoNgy  distributions  that  are  of  use  to  nearshore  modelers  seeking  boundary 
conditions  for  nearshore  wave  transformation  processes  or  the  results  of  near¬ 
shore  wave  generation  processes.  As  discuss^  in  Chapter  3,  rather  distinc¬ 
tive  spectral  shapes  arise  from  this  analysis  that  suggest  important  surf  zone 
processes  during  periods  of  high  wave  energy. 
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2  Data  Set 


Experiment  Site 

Data  used  in  this  study  were  acquired  at  the  Field  Research  Facility  (FRF) 
of  the  Coastal  Engineering  Research  Center  (CERC),  U.S.  Army  Engineer 
Waterways  Experiment  Station.  The  FRF  is  located  on  the  Outer  Banks  of 
North  Carolina  near  the  village  of  Duck  (Figure  1).  The  layout,  function,  and 
capabilities  of  the  FRF  are  described  by  Birkemeier  et  al.  (1985).  Of  rele¬ 
vance  to  this  study  are  the  continental  ^elf  bathymetry,  which  acts  to  steer 
the  waves  that  propagate  across  it,  and  the  general  nature  of  the  local  wind 
climate,  which  generates  the  higher-frequency  waves  observed  at  this  site. 

As  shown  in  Figure  1,  the  coastline  in  the  vicinity  of  the  FRF  is  nearly 
straight  for  several  tens  of  kilometers  north  and  south.  It  is  oriented  such  that 
a  shore-normal  line  (directed  seaward)  is  very  nearly  70  deg  from  true  north. 
Waves  and  onshore  winds  can  approach  this  site  along  an  easterly  ISOnleg  arc 
from  340  deg  to  160  deg  true,  llie  adjacent  continental  shelf  is  wide,  rela¬ 
tively  shallow,  and  of  somewhat  complex  bathymetry.  The  direction  of  near¬ 
est  ^proach  of  the  100-m  isobath  (approximately  the  SO-fathom  line  in  Figure 
1),  which  indicates  the  shelf  break,  is  tq)proximately  due  east,  or  about  20  deg 
south  of  the  shore-normal  line,  and  about  80  km  distant.  Long  wind  waves 
and  swell  generated  in  the  deep  Atlantic  Ocean  tend  to  refract  at  the  shelf 
break  and  :q)proach  the  FRF  site  from  this  direction.  A  typical  bottom  slope 
for  the  shelf  is  1  m/km,  but  this  is  interrupted  by  numerous  features  of  1-  to 
10-km  horizontal  scales  and  10-m  vertical  scales  positioned  irregularly  across 
the  shelf.  These  features  act  primarily  through  reifaction  as  wave  scattering 
lenses,  and  may  account  for  some  of  the  broad  directiotud  spreading  observed 
at  this  site  (Long  and  Oltman-Shay  1991). 

The  site  is  subject  to  a  variety  of  wind  conditions,  which  give  rise  to  a 
diverse  set  of  directional  wave  conditions.  Windless  days  are  rare.  Light 
winds  with  speeds  less  than  about  10  m/sec  are  very  common,  and  produce 
short,  high-frequency  waves  of  low  intensity,  usually  on  a  likewise  low-inten¬ 
sity  background  swell  arriving  from  the  deep  Atlantic.  For  the  data  set  used 
in  this  study,  the  primary  source  of  high-energy  waves  is  the  wind  associated 
with  passing  fronts,  many  of  which  have  embedded  low-pressure  cells.  Fron¬ 
tal  systems  are  usually  oriented  northeast-southwest  and  migrate  from  west  to 
east.  Frontal  winds  thus  tend  to  begin  by  blowing  from  the  southwest,  and 
then  shift  to  the  northeast  as  the  front  passes.  Though  the  southwest  winds 
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Figure  1 .  FRF  location  and  offshore  bathymetry 


have  an  offshore  blowing  component  that  tends  to  oppose  arriving  waves, 
there  is  also  a  southerly  component  that  generates  waves  propagating  north¬ 
ward  along  the  coast.  The  shallow  nearshore  bathymetry  tends  to  reft'act  these 
waves  onshore.  Consequently,  newly  generated  waves  from  the  southeast  are 
frequently  observed  when  a  frontal  system  approaches  the  FRF  site. 

Another  source  of  high-energy  waves  is  the  long-period,  well-sorted  swell 
that  radiates  from  distant  hurricanes.  Though  no  hurricanes  passed  over  the 
FRF  during  the  period  of  this  study,  several  made  landfall  well  south  of  this 
site  and  others  have  migrated  northward  past  this  site,  but  well  offshore. 
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These  waves  are  distinguished  from  those  associated  with  fronts  in  that  wave 
energy  is  distributed  over  narrow  ranges  of  frequency  and  direction. 


Directional  Wave  Gauge 

The  primary  instrument  used  in  this  study  is  a  high-resolution,  linear  array 
directional  wave  gauge.  It  consists  of  two  parts:  an  array  of  sensors  that 
sample  sea-surface  displacement  at  several  points  in  (horizontal)  space  and  an 
algorithm  that  converts  these  data  into  estimates  of  the  frequency  and  direc¬ 
tional  distribution  of  wave  energy.  The  array  consists  of  nine  pressure  gauges 
mounted  near  the  bottom  along  the  8-m  isobath  about  900  m  offshore  and  to 
the  north  of  the  FRF  research  pier  (Figure  2).  Spacing  between  the  gauges 
follows  the  array-design  guidance  of  Davis  and  Regier  (1977).  The  8-m  wat^ 
depth  is  shallow  enough  that  relatively  noise-ffee  wave  signals  can  be  detected 
up  to  a  frequency  of  about  0.32  Hz.  Waves  of  that  frequency  in  8  m  of  water 
have  a  wavelength  of  about  IS  m.  The  minimum  gauge  spacing  in  the  array 
is  S  m  and  spacings  increment  in  units  of  S  m,  so  that  these  high-frequency 
waves  are  well-sampled.  The  overall  length  of  the  array  is  255  m. 

Signals  from  each  of  the  pressure  sensors  are  sampled  at  2  Hz  for  2  hr 
16  min  (16,384-point  data  records).  Pressure  signals  are  windowed  in  IS 
half-l2q>ped  segments  of  17  min  duration  and  then  digitally  Fourier  trans¬ 
formed  to  the  frequency  domain.  Transformed  pr^sure  record  segments  are 
converted  to  transformed  sea-surface  displacement  record  segments  out  to  the 
0.32-Hz  cut-off  frequency  using  the  pressure  response  function  of  linear  wave 
theory.  Raw  cross-spectral  estimates  for  each  unique  pair  of  gauges  are  then 
computed  for  each  segment.  Smooth  cross-spectral  estimates  are  conq}uted  by 
averaging  over  the  15  segments  and  then  furAer  averaging  over  10  adjacent 
frequency  bands.  The  frnal  discrete  frequency  bandwidth  is  0.(X)977  Hz  and 
each  smoothed  cross-spectral  estimate  has  about  200  degrees  of  freedom  fol¬ 
lowing  Welch’s  method  as  described  by  Oppenheim  and  Schafer  (1975). 

Wave  energy  directional  distributions  are  found  by  forming  a  matrix  of 
cross-spectral  estimates  for  each  discrete  frequency  band  and  then  ^plying  the 
Iterative  Maximum  Likelihood  Estimation  (IMLE)  algorithm  derived  by 
Pawka  (1982,  1983)  as  a  refinement  to  the  ocean  wave  Maximum  Likelihood 
Estimation  method  described  by  Davis  and  Regier  (1977).  The  convergence 
parameters  used  by  Pawka  (1983)  were  used  here  identically.  There  is  no 
formal  estimate  of  error  for  IMLE,  and  tests  of  verity  must  be  done  on  a  case- 
by-case  basis  with  synthetic  data.  Preliminary  t^ts  of  the  FRF  array  indicate 
that  peaks  of  directional  distributions  can  be  resolved  to  within  about  5  deg 
and  that  spectral  estimates  within  a  directional  distribution  vary  with  a  stan¬ 
dard  deviation  of  about  30  percent  when  subjected  to  random  noise  appropri¬ 
ate  for  200  degrees  of  freedom.  Energy  estimates  for  this  study  were  comput¬ 
ed  in  discrete,  2-deg  arcs  for  the  full  180-deg  wave  approach  horizon. 

The  final  result  of  these  computations  is  a  regular  gridded  array  of  wave 
energy  density  as  a  function  of  discrete  frequency  and  direction.  The  frequen- 
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Figure  2.  FRF  linear  array  and  nearshore  bathymetry 

cy  dimension  is  in  28  bands  of  width  0.00977  Hz  covering  the  general  wind 
wave  frequency  range  of  0.05  to  0.32  Hz.  The  direction  dimension  is  in  91 
discrete  bins  of  2-deg  width  covering  the  180-deg  range  of  incident  wave 
approach  angles.  One  of  these  arrays  is  generated  each  time  a  data  acquisition 
process  known  as  a  collection  occurs. 
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Data  collection  at  the  FRF  occurs  at  least  four  times  daily,  at  0100,  0700, 
1300,  and  1900  Eastern  Standard  Time  (EST).  At  times  of  high  energy, 
typic^ly  when  characteristic  wave  heights  exceed  2  m,  additional  collections 
occur  at  0400,  1000,  1600,  and  2200  EST.  With  the  exception  of  the  first 
two  months,  Ais  collection  pattern  was  followed  for  the  full  duration  of  this 
riata  set,  extending  from  September  1986  through  August  1991.  September 
and  October  1986  were  the  months  of  the  SUPERDUCK  experiment 
(Crowson  et  al.  1988),  during  which  time  collections  were  3  hr  59  min  in 
duration,  centered  at  times  of  high  and  low  tide.  For  the  frequency-direction 
spectra,  these  long  records  were  processed  as  two  collections  of  2  hr  16  min, 
with  an  overlapping  34  min  of  data  shared  between  the  two  collections. 

Excepting  occasional  breaks  for  array  maintenance  and  repair,  and  two 
gaps  of  several  months,  data  coverage  was  continuous.  One  gap,  from  No¬ 
vember  1986  to  January  1987,  occurred  after  SUPERDUCK  and  is  omitted 
from  this  analysis  because  data  collections  were  reduced  to  34  min  in  dura¬ 
tion,  too  short  to  be  statistically  consistent  with  the  longer  records  used  in  the 
bulk  of  this  study.  The  second  gsq)  extended  from  December  1987  through 
August  1988,  beginning  when  a  fleet  of  fishing  trawlers  seriously  damaged  the 
array  and  ending  after  repair  parts  were  obtained  and  sufficient  fair  weather 
was  available  for  extensive  diving  operations. 

Altogether,  the  working  data  set  used  in  this  study  consists  of  6,759  dis¬ 
crete  frequency-direction  spectra.  Because  part  of  this  study  involves  group¬ 
ing  and  averaging  various  members  of  this  data  set,  it  is  noted  that  some 
aliasing  will  be  present  due  to  the  data  gaps  and  the  odd  sampling  pattern  dur¬ 
ing  SUPERDUCK.  Aliasing  in  this  context  means  that  the  time  axis  of  the 
sampling  is  not  uniform.  Consequently,  some  classes  of  spectra  may  be  over- 
or  under-represented.  This  condition  is  considered  acceptable  for  this  prelimi¬ 
nary  study  where  the  intent  is  to  employ  the  maximum  possible  number  of 
observations.  Additional  aliasing  by  higher  sampling  rates  during  high  wave 
conditions  is  not  a  problem  in  this  study  because  wave  height  is  one  of  the 
grouping  parameters  and  the  2-m  criterion  for  higher  sampling  occurs  at  one 
of  the  height  class  boundaries. 


Parameter  Definitions 

The  main  objective  of  this  study  is  to  classify  and  group  the  observed 
spectra  by  characteristic  wave  height,  period,  and  direction.  Many  ways  to 
define  rq>resentative  parameters  with  these  names  are  found  in  the  literature, 
especially  where  frequency-direction  spectra  are  used.  The  specific  parame¬ 
ters  used  in  this  study  are  defined  here. 

Figure  3  is  representative  of  spectra  obtained  from  the  FRF  array.  The 
three-dimensional  surface  depicts  the  discrete  frequency-direction  spectrum 
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FRF  Frequency-Direction  Spectrum 
for  2200  EST  on  1  Mar  91 
Hmo  “  0.80  m  fp  “  0.16  Hz  9p  "  -44  deg 


Figure  3.  Representative  frequency-direction  spectrum 


denoted*  ,  where  is  the  n*  of  -  28  discrete  frequencies  (f)  diat 
increment  in  units  of  4^  =  0.00977  Hz,  and  6^  is  die  m*  of  ilf  =  91  dis¬ 
crete  directions  that  increment  in  units  of  A0  =  2  deg.  Direction  is  defined 
such  that  6  represents  the  direction  from  which  waves  arrive,  measured  coun¬ 
terclockwise  from  shore-normal.  The  frequency  spectrum  in  the  right  panel  of 
Figure  3  is  the  conventional  spectrum  that  one  would  obtain  widi  data  from  a 
single  wave  staff  or  surface-corrected  pressure  gauge.  The  frequency  spec¬ 
trum  is  related  to  the  frequency-direction  spectrum  by 


if 


SK)  -  E  svjj  A» 

•••I 


(1) 


'  For  convenience,  lymboli  and  abbreviations  are  Usted  in  the  notation  (Appendix  D). 
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The  direction  spectrum  in  the  left  panel  of  Figure  3  is  the  directional  analogue 
of  the  frequency  spectrum,  representing  the  total  energy  in  each  direction  bin. 
The  direction  spectrum  is  related  to  the  frequency-direction  spectrum  by 


s(ej  .  £  at 

Because  the  frequency-direction  spectra  are  to  be  classified,  characteristic 
parameters  representing  wave  height,  frequency,  and  direction  are  deduced 
from  the  thre^imensional  part  of  Figure  3. 


Wave  height 

In  this  study,  characteristic  wave  height  is  the  conventional  energy-  or 
spectrum-based  height  .  It  is  defined  as 

/f  =  4<y  (3) 

MO 


where  a  is  the  standard  deviation  of  sea  surface  displacement,  which  is  tiie 
square  root  of  the  variance  of  sea  surface  displacement  within  the  bounding 
frequency  and  direction  limits  employed  in  spectral  computation.  In  terms  of 
the  frequency-direction  spectrum,  the  variance  is 

»  E  E  S(fJJ  AO  Af  (4) 

nal  M«1 


Frequency 

Charaaeristic  frequency  is  defined  here  as  the  frequency  associated  with 
the  peak  or  maximum  value  of  .  A  more  conventional  definition  is  the 
frequency  associated  with  the  peak  of  S(f^,  primarily  because  there  are  many 
more  observations  of  frequency  spectra  than  of  frequency-direction  spectra. 
These  two  peak  frequencies  can  differ  if  wave  energy  is  concentrated  near 
some  direction  at  one  frequency  and  more  directionally  spread  at  a  second 
frequency.  If  the  frequency  spectrum  shows  equal  amounts  of  energy  at  these 
two  frequencies,  the  frequency-direction  spectrum  will  show  a  strongs  peak 
at  the  first  frequency,  and  will  be  classified  with  equal  to  the  first 

frequency. 

This  distinction  is  important  because  the  more  conventional  definition  of/^ 
would  associate  the  frequency-direction  spectrum  with  the  second  frequency, 
and  would  subsequently  alter  statistical  properties  of  groups  of  spectra  identi¬ 
fied  by  classes  of  frequency.  The  definition  of  ^  used  in  this  study  is  justi- 
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fied,  however,  because  similarities  in  shapes  of  die  three-dimensional  spectra 
are  sought  here.  The  frequency  for  which  energy  is  most  concentrated  in  the 
three-dimensional  spectrum  is  Aus  the  appropriate  param^o-.  It  is  also  most 
appropriate  in  engineering  studies  where  directionality  and  groupiness  proper¬ 
ties  are  included,  because  waves  of  extreme  height  are  more  likely  to  occur  at 
this  frequency. 


Direction 

In  the  same  manner  and  for  the  same  reasons  as  characteristic  frequency  is 
defined,  characteristic  direction  0^  is  identified  as  the  direction  at  which 
is  a  maximum.  It  is  noted  that  in  some  numerical  models,  like  CERC’s  Wave 
Information  Studies  (Corson,  Resio  and  Vincent  1980),  and  in  the  treatment  of 
data  from  some  low-resolution  directional  wave  gauges  (Howell  1992),  a  char¬ 
acteristic  direction  is  identified  as  the  mean  direction,  which  is  the  zeroth 
moment  of  a  directional  distribution  function,  usually  that  occurring  at  the 
frequency  of  the  peak  of  S(Q .  This  second  definition  can  vary,  sometimes 
dramatically,  from  the  6^  used  in  the  present  study,  but,  again,  the  present 
definition  is  simpler  and  assures  that  constituent  spectra  in  a  direction  class  all 
have  the  greatest  concentration  of  energy  at  the  same  direction.  It  is  likely 
that  future  studies,  involving  model  testing  or  gauge  intercomparisons,  will 
consider  mean  direction  as  a  characterizing  parameter. 


Classification  Scheme 


With  the  characteristic  parameters  defined,  the  complete  set  of  frequency- 
direction  spectra  can  be  classified  in  a  parameter  space  having  axes  , 

and  statistical  properties  of  constituent  spectra  within  a  class  can  be  computed. 
To  gain  statistical  confidence,  and  to  maintain  a  reasonable  overall  numbw  of 
resulting  classes,  tiie  full  ranges  of  each  of  the  three  parameters  are  subdivid¬ 
ed  into  small  numbers  of  intervals  so  that  any  spectrum  can  be  identified  with 
a  specific  height,  frequency,  and  direction  class. 


Height  classes 

Within  the  present  data  set,  varied  from  0.15  m  to  4.70  m.  Two 
constraints  affect  how  wave  hei^t  class  boundaries  are  defined.  First,  it  is 
wave  energy  that  a  spectrum  represents,  so  that  classification  should  be  by 

energy  content,  which  varies  as  and  not  just  wave  height.  Second,  there 
are  far  more  observations  at  small  than  at  large,  suggesting  that  class 
boundaries  should  be  closer  together  at  small  wave  heights  to  get  a  more  even 
number  of  cases  in  each  of  the  height  classes.  Both  of  these  constraints  are 
satisfied  by  using  a  geometric  distribution  of  class  boundaries  instead  of  a 
simple  linear  distribution. 
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A  geometric  distribution  means  that  the  ratio  of  one  height  class  boundary 
to  the  previous  height  class  boundary  is  constant.  When  squared  to  rq)resent 
energy,  this  scheme  results  in  a  constant  percentage  change  of  energy  across 
all  classes.  In  this  study,  eight  height  classes  were  used.  The  geometric 
distribution  results  in  nine  height  class  boundaries,  at  «  O.IS,  0.23,  0.35, 
0.55,  0.84,  1.29,  1.99,  3.06,  and  4.70  m.  As  noted  before,  one  of  the  boun¬ 
daries  occurs  very  near  to  2  m  so  that  some  of  the  biasing  effects  of  higher 
spectral  sampling  ^ve  this  wave  height  are  isolated  to  the  two  highest  wave 
height  classes. 


Frequency  desses 

All  spectra  were  defined  in  terms  of  28  discrde  frequency  bands  ranging 
from  about  0.054  Hz  to  about  0.318  Hz.  Hiis  range  is  divided  into  seven 
classes,  each  having  a  nominal  bandwidth  of  about  0.04  Hz.  All  classes 
contain  four  possible  discrete  peak  frequencies  except  the  lowest-frequency 
class,  which  contains  only  three.  The  class  boundaries  occur  at  /  -  0.054, 
0.083,  0.122,  0.161,  0.200,  0.240,  0.279,  and  0.318  Hz.  For  plotting  pur¬ 
poses,  the  classes  are  identified  by  nominal  band-center  frequencies  of  0.06, 
0.10,  0.14,  0.18,  0.22,  0.26,  and  0.30  Hz. 


Direction  classes 

The  180-deg  arc  that  represents  the  possible  range  of  is  divided  into  19 
bins.  Each  bin,  except  the  first  and  last,  has  a  range  of  10  deg,  and  is  cen¬ 
tered  on  an  azimuth  that  is  an  integer  multiple  of  10  deg  from  shore-normal. 
The  first  and  last  bins  only  have  ranges  of  5  deg  because  they  terminate  at 
90  deg  and  -90  deg,  respectively.  The  bin  boundaries  thus  follow  tihe  pattern 
90  deg,  85  deg,  75  deg,  65  deg,...,  5  deg,  -5  deg,...,  -65  deg,  -75  deg, 

-85  deg,  -90  deg.  For  plotting  purposes,  direction  classes  are  assigned  nomi¬ 
nal  values  of  90,  80,  70,...,  10,  0,  -10,...,  -70,  -80,  and  -90  deg.  Though  die 
first  and  last  bins  are  only  5  deg  wide,  it  turns  out  that  no  spectra  fall  into 
those  classes,  so  the  analysis  is  not  compromised. 


Class  Distribution  of  Spectra 

With  eight  wave  height  ranges,  seven  frequency  ranges,  and  nineteen 
direction  ranges,  there  are  a  total  of  1,064  possible  classes  with  which  to 
characterize  the  6,759-case  spectral  database.  If  these  spectra  were  uniformly 
distributed  across  all  possible  classes,  there  would  be  roughly  six  cases  per 
class,  enough  to  do  some  sinqile  statistics.  Of  course,  the  spectra  are  not 
distributed  uniformly,  so  it  is  of  interest  to  see  how  the  spectra  are  distributed. 
If  the  total  number  of  spectra  that  fall  into  a  given  class  is  divided  by  the  total 
number  of  cases  in  the  database  and  multiplied  by  100  percent,  a  set  of  empir¬ 
ical  distribution  functions  can  be  computed. 
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Figure  4  is  one  way  to  display  the  three-dimensional  distribution  of  cases 
in  the  wave  height,  peak  frequency,  and  peak  direction  domains.  There  is  one 
subplot  for  each  wave  height  class.  Within  each  subplot  is  a  diree-dimen- 
sional  surface  that  shows  as  a  percentage  the  number  of  all  cases  that  fall  into 
a  particular  frequency  and  direction  class.  The  header  for  each  subplot  shows 
the  percentage  of  all  cases  that  fall  into  the  corresponding  wave  height  class. 
Side  panels  in  each  subplot  show  percentages  of  all  cases  associated  with  a 
particular  frequency  class  (right  panels)  or  direction  class  Oeh  panels)  within  a 
wave  height  class.  All  parts  of  ^1  subplots  are  drawn  to  the  same  scale. 

The  wave  height  class  with  the  lowest  energy  has  too  few  cases  to  be 
visible  at  the  scale  shown,  but  also  represents  so  little  energy  as  not  to  be 
important  dynamically.  The  next  three  wave  height  classes  show  two  impor¬ 
tant  features  of  the  wave  climate  at  the  FRF.  Chie  is  the  pervasive,  low-  to 
mid-frequency  swell  that  tends  to  arrive  from  south  of  shore-normal  (negative 
peak  directions),  suggesting  generation  somewhere  in  the  deep  Atlantic  and 
refraction  at  the  shelf  break  (Figure  1).  The  pervasiveness  is  reflected  in  the 
high  percentage  of  cases  represented  by  these  classes.  The  other  feature  is  the 
set  of  spectra  associated  with  the  initial  stages  of  wind  wave  generation. 

These  cases  tend  to  have  peaks  at  middle  to  high  frequencies  and  are  distin¬ 
guished  by  having  high  approach  angles  from  both  the  north  (initial  stages  of 
northeasterly  storms)  and  the  south  Ofossibly  refracted  waves  generated  by  the 
southwesterly  winds  of  approaching  fronts).  Spectra  with  high  frequency 
peaks  very  seldom  occur  at  near-normal  incidence. 

At  still  higher  energies,  there  are  three  dominant  characteristics.  First, 
there  are  very  few  high-frequency  peaks,  indicating  continued  sea-state 
growth.  Second,  there  are  fewer  peak  directions  in  the  southerly  quarter, 
showing  the  increasing  dominance  of  northeast  winds  in  wave  generation. 
Third,  the  most  frequent  cases  are  of  low  peak  frequency  and  more  centered 
at  shore-normal  incidence,  indicating  the  presence  of  refracted,  locally  gener¬ 
ated  long  waves.  The  two  classes  containing  the  hipest  wave  heights  do  not 
contain  a  high  percentage  of  cases,  but  are  very  in^rtant  dynamically  be¬ 
cause  they  represent  the  peak  conditions  of  storms  diat  incur  great  changes  to 
nearshore  bafoymetry,  increased  beach  erosion,  and  structural  damage.  These 
two  classes  are  characterized  by  low-frequency  peaks  at  normal  to  slightly 
north  of  normal  incidence. 

This  simple  interpretation  of  the  distribution  surfaces  of  Figure  4  is  consis¬ 
tent  with  identifiable  physical  processes,  and  suggests  that  meaningful  insight 
can  be  gained  by  a  three-parameter  characterization  of  frequency-direction 
spectra.  However,  it  is  6a  from  complete.  Considerable  variation  in  the 
structure  of  spectral  shs^es  within  a  wave  height,  peak  frequency,  and  peak 
direction  classification  is  quite  probable.  These  variations  can  be  very  impor¬ 
tant  dynamically,  so  the  next  step  is  to  examine  frequency-direction  spectral 
sh^e  within  individual  three-parameter  classifications. 
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Figure  4.  Distribution  of  spectral  cases  in  parameter  space 
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3  Characteristic  Spectra 


More  insight  into  spectral  structure  beyond  the  gross  climatological  inter¬ 
pretation  afforded  by  simply  counting  cases  is  gained  by  examination  of  statis¬ 
tical  properties  of  constituent  spectra  within  each  wave  height,  peak  frequen¬ 
cy,  and  peak  direction  class.  Of  particular  interest  are  the  frequency  and 
directional  energy  distributions  that  are  typical  or  characteristic  of  each  classi¬ 
fication.  The  goal  is  to  provide  meaning^  definitions  of  characteristic  spec¬ 
tra  as  whole  entities,  thus  extending  the  work  of  Long  and  Oltman-Shay 
(1991),  who  examined  only  the  first  year  of  data,  and  who  looked  only  at 
directional  distributions  of  energy  independently  of  the  frequency-direction 
spectral  context  from  which  they  came. 


Mean  Spectra 

Perhaps  the  easiest  way  to  obtain  a  characteristic  spectrum  from  a  s^  of 
constituent  spectra  is  to  compute  an  average  spectral  density  for  each  frequen¬ 
cy  band  and  direction  bin  in  the  frequency-direction  plane.  If  there  are  K 
constituent  spectra  in  a  particular  height-period-direction  class,  a  first  attenq>t 

at  an  average  spectrum  can  be  defined  mathematically  as 

-  4  E  ® 

whne  is  the  k*  constituent  spectrum  in  the  class. 

If  a  mean  spectrum  can  be  computed,  it  is  reasonable  that  some  measure  of 
variability,  such  as  the  standard  deviation,  can  also  be  computed.  Variability 
among  constituent  members  of  a  class  of  spectra  can  arise  from  two  primary 
causes.  One  cause  is  variation  in  spectral  sbtq>e  from  such  properties  as  sec- 
oiKlary  spectral  peaks  or  the  amount  of  directional  spread.  The  other  cause  is 
variation  in  energy  content.  Even  though  constituent  spectra  are  classifred  by 
energy  content  by  way  of  H^,  tiiere  is  a  finite  range  of  energies  within  a 
class.  Because  of  the  geometric  spacing  of  height  class  boundaries,  the  total 

energy,  which  is  proportional  to  can  vary  by  a  factor  of  about  2.4  across 
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one  height  class.  The  present  interest  is  in  variability  of  spectral  shiy)e,  so  it 
is  desired  to  minimize  the  influence  of  energy  content  in  the  analysis.  One 
way  to  isolate  spectral  shape  from  energy  content  is  to  normalize  each  constit¬ 
uent  spectrum  by  its  own  volume,  which,  by  Equation  2,  is  the  variance  of 
sea  surface  displacement  o^.  Normalized  spectra  are  defined  by 

(6) 

Ok 


where  a\  is  the  variance  contained  in  the  k"*  constituent  spectrum  in  a  class. 
Normalized  spectra  are  dimensionless  and  have  unit  v  olume,  so  that  compar¬ 
ing  normalized  spectra  emphasizes  spectral  shape,  as  desired.  Note  that  by 

Equations  3  and  4,  al  •  where  is  the  characteristic  wave 

height  associated  with  the  constituent  spectrum  . 

The  remainder  of  the  analysis  discussed  in  this  rqK)rt  will  be  based  on 
normalized  spectra.  Though  the  normalized  spectra  are  individually  indepen¬ 
dent  of  energy  content,  classification  by  characteristic  wave  height  is  retained 
because,  as  will  be  seen,  spectral  shapes  associated  with  low  energy  are  dis¬ 
tinctively  different  from  those  with  high  energy. 

The  mean  of  a  set  of  K  constituent  normalized  spectra  within  a  class  is 

-  i  t  sitfM  f 

A  Jt«t 


Standard  Deviation  Spectra 

To  show  the  variability  in  spectral  sht^e,  a  standard  deviation  is  conq>uted 
for  each  frequency  band  and  direction  bin  of  the  constituent  spectra.  From  K 

normalized  spectra,  the  result  is  a  standard  deviation  spectrum  com¬ 

puted  following  the  conventional  definition  of  standard  deviation  by 


si(fM 


t 

K  -  1 


(8) 


Computation  of  a  standard  deviation  requires  a  set  of  more  than  one  constitu- 
rat  spectrum.  In  this  study,  mean  and  standard  deviation  spectra  are  comput¬ 
ed  o^y  for  those  classes  having  at  least  three  constituent  spectra.  Con¬ 
sequently,  some  of  the  very  sparsely  populated  classes  are  eliminated. 
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The  standard  deviation  spectrum  is  a  surface  that  shows  variability  of 
normalized  spectral  densities  as  a  function  of  frequency  and  direction.  A 
useful  parameter  that  gives  a  bulk  measure  of  spectral  variability  is  the  volume 
of  the  standard  deviation  spectrum  defined  by 

N  M 

‘'.-EE  ^  AT 

n^l  mal 


When  is  small,  the  volume  of  the  standard  deviation  spectrum  is  small, 
and  the  interpretation  is  that  the  constituent  spectra  making  up  the  mean  spec¬ 
trum  are  all  rather  alike.  When  V  is  large,  the  volume  of  the  standard  devia¬ 
tion  spectrum  is  large,  suggesting  ^at  the  constituent  spectra  have  rather 
different  shapes. 


Sample  Case 

Figure  S  illustrates  the  mean  (upper  graph)  and  standard  deviation  (lower 
graph)  spectra  for  one  of  the  height-frequency-direction  classes,  the  bounding 
parameters  for  which  are  shown  in  the  plot  header.  Representative  frequency 
and  direction  spectra,  computed  using  Equations  1  and  2,  respectively,  are 
shown  in  the  side  panels  of  the  plots. 

The  mean  spectrum  in  Figure  S  is  characterized  by  a  broad,  rather  skewed 
directional  distribution  near  the  peak  frequency,  and  narrow,  relatively  low- 
energy  directional  distributions  at  the  higher  frequencies.  Tiis  shape,  which 
is  the  mean  or  expeaed  shape  from  95  constituent  spectra,  is  only  grossly 
^proximated  by  a  unidirectional,  monochromatic  wave  field  having  the  same 
peak  frequency  and  peak  direction.  A  unidirectional,  monodiromatic  wave 
field  would  be  represented  in  spectral  form  in  Figure  S  as  a  spike  having  unit 
volume  with  base  dimensions  A/  x  Ad  centered  at  the  peak  frequency  and 
peak  direction  of  the  mean  spectrum.  Lacking  other  knowledge  beyond  the 
height,  period,  and  direction  characterizing  parameters,  a  unidirectional, 
monochromatic  spectrum  would  be  the  only  way  to  represent  this  class  of 
spectra.  As  will  be  seen  in  Chapters  4  and  5,  direction  sensitive  wave  pro¬ 
cesses,  such  as  energy  flux  and  radiation  stress,  tend  to  have  very  differrat 
values  when  comput^  from  a  unidirectional,  monochromatic  spectrum  as 
compared  to  a  representative  mean  spectrum  like  that  shown  in  Figure  S. 

Of  some  concern  is  the  size  and  shape  of  the  standard  deviation  spectrum 
shown  in  Figure  5.  Both  the  mean  and  standard  deviation  spectral  gnqihs  are 
drawn  to  the  same  scales,  and  it  is  evident  that  the  standard  deviation  spec¬ 
trum  is  shaped  somewhat  like  the  mean  spectrum  and  is  about  the  same  size, 
as  can  be  seen  visually  and  as  measured  by  V^,  which  is  slightly  largo:  than 
unity.  It  would  be  preferred  that  the  standard  deviation  spectrum  be  small  at 
all  frequencies  and  directions,  indicating  that  the  mean  spectrum  is  clearly 
r^resentative  of  its  class,  but  this  evidently  is  not  the  case  here. 
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Mean  and  Standard  Deviation  of  Normalized  Spectra 
95  Cases  0.55  <  (m)  <  0.84 

0.08  <  fp  (Hz)  <  0.12  -35  <  6^  (deg)  <  -25 

Mean  Volume  «■  1.01  Standard  Deviation  Volume  —  1.10 


Figure  5.  Examples  of  mean  and  standard  deviation  spectra 
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There  are  two  probable  main  reasons  why  the  standard  deviation  spectrum 
is  so  large.  One  reason  is  that  there  may  be  constituent  spectra  having  the 
general  shape  of  the  mean  spectrum  in  Figure  S.  but  displaced  slightly  along 
the  frequency  or  direction  axes  within  the  bounding  frequencies  and  directions 
of  the  classification  scheme.  In  the  classification  used  here,  each  direction 
class  encompasses  five  of  the  discrete  direction  bins  that  define  a  constituent 
spectrum,  a^  the  peak  direction  of  a  constituent  spectrum  could  fall  in  any 
one  of  those  five  bins.  Likewise,  the  frequency  classes  include  four  discrete 
frequency  bands  that  define  a  constituent  spectrum,  and  the  peak  frequency  of 
a  constituent  spectrum  could  reside  in  any  one  of  those  bands.  Collectivdly, 
within  a  class,  there  are  20  possible  unique  frequency-direction  locations  that 
could  represent  a  constituent  spectral  pe^.  When  the  constituent  spectra  are 
steq>-sided,  the  standard  deviation  computation  could  involve  a  number  of 
small  values  and  a  number  of  large  values,  depending  on  how  the  constituent 
spectral  peaks  are  distributed  in  the  class  ranges  of  frequency  and  direction, 
mid  this  can  give  rise  to  a  large  standard  deviation.  When  this  reason  is  the 
cause  of  a  large  standard  deviation,  the  resulting  mean  spectrum  is  a  quite 
reasonable  representation  of  the  constituent  members,  being  somewhat 
smoother  and  of  slightly  broader  spread,  but  still  rmaining  the  fundamental 
constituent  shape. 

The  other  reason  is  that  the  constituent  spectra  may  have  significant  struc¬ 
tural  differences  that  are  not  reflected  by  the  characteristic  frequency  and 
direction  designations.  These  differences  may  include  broad  versus  narrow 
directional  or  frequency  distributions  of  energy,  or  die  existence  of  significant 
secondary  peaks  or  side  lobes  in  some  of  the  constituent  spectra  that  are  lack¬ 
ing  in  others.  This  type  of  behavior  can  be  severe  in  that  the  mean  spectrum 
does  not  represent  any  one  of  the  constituent  sp^rtra  very  well. 

That  both  reasons  are  extant  in  the  present  analysis  is  made  evident  by 
Figure  6,  which  shows,  in  somewhat  compact  form,  12  of  the  95  constituent 
spectra  that  make  up  the  mean  and  standard  deviation  spectra  of  Figure  5. 
Most  of  the  individual  spectra  in  Figure  6  have  broad,  somewhat  skewed 
directional  distributions  near  the  low-frequency  peaks  and  at  least  the  begin¬ 
nings  of  directionally  nanow,  high-frequency  tails,  like  the  mean  spectrum  of 
Figure  S.  However,  the  top,  left-most  spectrum  clearly  has  secondary  struc¬ 
ture  in  the  forms  of  a  secondary  peak  at  a  frequency  a^  a  direction  well  away 
from  the  primary  peak  used  to  classify  this  case,  and  a  low-energy,  high- 
frequency  tail  along  a  set  of  directions  that  are  very  different  from  those  of 
the  high-frequency  tail  predominant  in  the  mean  spectrum  of  Figure  5.  Clear¬ 
ly,  the  mean  spectrum  of  Figure  5  does  not  represent  this  case  as  well  as  it 
does  the  oftier  cases  shown  in  Figure  6. 

One  way  to  deal  with  this  problem  is  to  make  further  distinctions  about 
details  of  spectral  structure,  essentially  by  adding  classifications  for  such 
things  as  characteristic  directional  spread  and  number  of  dominant  modes,  and 
redoing  the  analysis.  Such  a  procedure  is  beyond  the  scope  of  the  present 
enterprise,  however,  because  the  effort  here  is  to  ascertain  some  of  the  conse¬ 
quences  of  a  three-parameter  characterization  of  .frequency-direction  spectra. 
The  present  problem  is  one  such  consequence;  averages  of  constituent  spectra 
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monocbromatic  spectra.  A  quantification  of  this  statement  is  conq>uted  for 
some  specific  applications  in  Giapters  4  and  S. 


General  Set  of  Classes 

The  behavioral  properties  described  for  the  single  class  of  spectra  rq)re- 
sented  by  Figures  5  and  6  are  common  to  all  of  the  spectral  classes  analyzed 
in  dtis  study.  Mean  spectra  for  all  classes  having  at  least  three  constituent 
spectra  are  illustrated  in  Appendix  A.  Corresponding  standard  deviation 
spectra  are  illustrated  in  Appendix  B.  The  variability  of  constituent  spectra  as 
expressed  by  the  volumes  of  the  standard  deviation  spectra  V,  is  illustrated  for 
all  classes  in  Figure  7.  Empty  classes,  those  having  fewer  than  three  constitu¬ 
ents,  are  plotted  as  having  V,  =  0  in  Figure  7. 

For  classes  that  are  not  empty.  Figure  7  shows  that  the  V,  are  all  of  order 
one,  ranging  from  about  0.3  to  about  1.1,  and  tend  to  be  largest  for  spectra 
with  low  peak  frequencies  in  intermediate  wave  height  classes.  The  examples 
of  mean,  standard  deviation,  and  constituent  spectra  shown  in  Figures  5  a^  6 
represent  a  case  where  V,  is  large.  Of  inter^t  is  the  character  of  a  class 
where  V,  is  smaller. 


High-Energy  Case 

Figure  8  illustrates  the  mean  and  standard  deviation  spectra  for  a  class 
where  V,  is  about  half  that  for  the  class  illustrated  in  Figure  S.  The  class 
shown  in  Figure  8  is  one  of  the  most  important  classes  because  it  rq)resents 
conditions  of  very  high  energy.  The  mean  spectrum  is  characterized  by  a 
directionally  bro^,  low-frequency  peak  and  a  distinctive,  directionally  bimod- 
al,  high-frequency  tail.  This  general  structure  was  also  observed  by  Long  (in 
pr^aration)  in  a  study  of  the  evolution  of  directional  spectra  during  individual 
storm  events.  Though  derived  from  the  same  data  set,  the  present  result  is 
notable  because  it  suggests  that,  with  minor  variations,  the  same  general  struc¬ 
ture  appears  in  the  same  ranges  of  frequency  and  direction  for  many  storms. 

Spectral  variation  for  this  case  is  depicted  in  the  lower  graph  in  Figure  8. 
The  shiq>e  of  the  standard  deviation  spectrum  is  roughly  like  Aat  of  the  mean 
spectrum,  with  the  exception  of  one  strong  peak  at  low  frequency  and  cen¬ 
tered  near  the  shore-normal  direction.  The  source  of  this  peak  is  shown  in 
Figure  9,  which  shows  all  10  of  the  constituent  spectra  in  this  class.  The 
strong  peak  in  Figure  8  is  due  primarily  to  the  last  (right  spectrum,  bottom 
row)  constituent  spectrum  shown  in  Figure  9.  All  the  other  spectra  in  Fig¬ 
ure  9  are  of  quite  .similar  structure,  having  the  directionally  broad,  low- 
frequency  peaks  and  directionally  bimodal,  high-frequency  tails  of  the  mean 
spectrum  of  Figure  8.  Consequently,  most  of  the  variance  shown  in  the  lower 
graph  of  Figure  8  is  due  to  relatively  minor  variations  in  amplitude  and 
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Figure  8.  Mean  and  standard  deviation  spectra  for  a  high-energy  case 
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Discussion 


One  of  the  primary  results  of  this  study  is  the  s^  of  characteristic  spectra 
discussed  in  this  chi4)ter  and  illustrated  in  Appendix  A.  Analysis  of  variabili¬ 
ty  and  its  causes  suggests  that  the  characteristic  spectra  reflect  at  least  the 
gross  structures  of  the  constituent  spectra  in  a  class.  Hence,  the  charaaeristic 
spectra  are  much  to  be  preferred  over  unidirectional,  monochromatic  spectra 
for  engineering  computations  related  to  this  and  similar  sites  because  they  are 
more  like  the  spectra  of  naturally  occurring  wave  fields. 

A  modeler  seeking  seaward  boundary  conditions  for  a  nearshore  processes 
model  or  landward  results  for  a  model  of  wave  transformation  across  the 
continental  shelf  can  use  these  characteristic  spectra  as  guides  or  templates  of 
real  wave  properties.  A  model  that  produces  results  near  to  these  characteris¬ 
tic  spectra,  within  the  variability  quantified  in  this  chapter  and  in  the  standard 
deviation  spectra  illustrated  in  Appendix  B,  can  be  judged  as  being  valid 
within  the  range  of  natural  variability.  Consequently,  the  characteristic  mean 
and  standard  deviation  spectra  found  here  should  be  of  considerable  aid  in 
model  testing  and  development. 

Up  to  this  point  in  this  report,  it  has  only  been  emphasized  that  the  charac¬ 
teristic  spectra  found  here  are  significantly  different  in  shape  from  the  unidi¬ 
rectional,  monochromatic  wave  spectra  that  are  the  commonly  used  models  for 
a  three-parameter  defoiition  of  wave  conditions.  Quantitative  differences  in 
properties  used  in  modeling  and  engineering  computations  are  needed  to  em¬ 
phasize  the  importance  of  using  the  characteristic  spectra  instead  of  unidirec¬ 
tional,  monochromatic  waves,  or,  alternatively,  the  extent  to  which  estimates 
using  unidirectional,  monochromatic  waves  can  be  corrected  to  reflect  natural 
conditions  with  higher  fidelity.  Examples  of  these  types  of  computations  are 
given  in  the  following  two  chapters. 
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4  Longshore  Energy  Flux 


As  a  useful  example  of  the  sensitivity  of  engineering  conq>utations  to 
choice  of  spectral  definition,  longshore  energy  fluxes  are  estimated  and  com¬ 
pared.  The  relevance  of  longshore  energy  flux  is  in  its  use  as  part  of  a  model 
for  longshore  sediment  transport  contained  in  the  SPM  (1984).  That  model 
asserts  that  longshore  sediment  transport  is  directly  proportional  to  longshore 
energy  flux,  so  that  any  errors  in  estimating  longshore  energy  flux  translate 
directly  into  estimation  errors  of  longshore  sediment  transport. 


Definitions 

In  terms  of  a  discrete  frequency-direction  spectrum  ,  a  general 

expression  for  longshore  energy  flux  P^  is 

E  E  S(fJJ  cos  0,  sin  A0  Af 

M*!  m*! 


where  p  is  water  density,  g  is  gravitational  acceleration,  and  is  wave 

group  velocity  :q)propriate  to  frequency  /,  and  water  depth  d.  Equation  10 
expresses  the  amount  of  longshore-direaed  wind  wave  energy  flux  through  the 
point  of  spectral  observation  per  unit  of  longshore  distance. 

To  concentrate  on  the  directional  sensitivity  of  estimates  of  longshore 
energy  flux,  it  is  convenient  to  normalize  both  sides  of  Equation  10  with  the 

expression  pg<Py/gT.  Division  by  pg,  whidt  varies  insignificantly,  elimi¬ 
nates  the  leading  coefficient  of  Equation  10.  Division  by  which  varies 
only  nominally  due  to  a  tide  range  of  about  1  m  in  the  8-m  water  depth, 
results  in  a  dimensionless  group  velocity.  Division  by  converts  the  fre¬ 
quency-direction  spectrum  to  the  normalized  form  S' (f ^,9 J  defined  by  Equa¬ 
tion  6  and  used  in  Qiapter  3.  The  result  of  the  normalization  is  a  dimension¬ 
less  longshore  energy  flux  P/ ,  which  has  the  form 
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(11) 


"  "  C(f  d\ 

P'l  -EY:  cos  e.  sin  A0  at 

aal  *•1  ^gd 

In  this  chiqpter,  estimates  of  pI  are  made  from  three  spectral  rq)resenta- 
tions.  One  is  a  monochromatic,  unidirectional  r^resentation  for  which  enti¬ 
ties  are  identified  by  the  subscript  mu .  Aimther  is  the  set  of  estimates  fi'om 
observed  frequency-direction  spectra,  which  served  as  the  constituent  spectra 
in  the  definitions  of  characteristic  spectra  in  (Zh^ter  3.  Entities  derived  from 
observed  frequency-direction  spectra  are  identified  by  the  subscript  fd.  The 
third  stt  of  estimates  is  derived  from  the  characteristic  spectra  found  in  Chap¬ 
ter  3  and  displayed  in  Appendix  A.  Entities  associated  with  characteristic 
spectra  are  identified  with  an  overbar  on  the  subscript  fd. 

When  a  monochromatic,  unidirectional  spectrum  is  employed.  Equation  11 
is  reduced  to  a  much  simpler  form.  A  non^ized  monochromatic,  unidirec¬ 
tional  spectrum  in  discrete  frequency-dimension  space  has  a  non-zero  magni¬ 
tude  equal  to  VAdAf  at  only  one  discrete  frequency  f  ^  and  one  discrete 
direction  6^  >  6^.  Consequently,  there  is  only  one  non-zero  value  within  the 
double  summation  of  Equation  11,  and  the  normalized  monochromatic,  unidi¬ 
rectional  measure  of  longshore  energy  flux  is 


cos  0.  sin  6, 

F  P 


(12) 


When  both  sides  of  Equation  12  are  multiplied  by  pg<P^f^,  and  it  is  noted 

f  1  * 

that  0^  >  1 1  —  I  for  monochromatic,  unidirectional  waves,  where  H  is 

H2I 

actual  wave  nei^t,  the  dimensional  form  of  longshore  energy  flux  defined  in 
the  SPM  is  recovered,  i.e., 

I  cos  sin  (13 


Within  a  height-frequency-direction  classification,  there  are  K  occurrences 
of  the  parameter  set  .  For  the  k*  occurrence,  the  monodiromatic, 

unidirectional  estimate  of  longshore  energy  flux  is 

cos  9^  sin  9^  (14) 

M 
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where  $^j,,  and  are,  respectively,  the  peak  frequency,  peak  direction, 
and  water  d^th  associated  with  the  k*  case.  A  parameter  that  represents  the 
set  of  K  monochromatic,  unidirectional  estimates  within  a  class  as  a  whole  is 
the  mean  of  the  set,  defined  by 


(15) 


A  standard  deviation  of  the  set  of  monochromatic,  unidirectional  estimates 
within  a  class  could  be  defined  here,  but  it  has  been  found  to  be  small  relative 
to  other  measures  of  variability.  Consequently,  it  is  not  used  in  this  analysis. 

When  dealing  with  a  group  of  constituent  frequency-direction  spectra  with- 

“  *  classification,  the  k*  spectrum  produces  an  estimate  of  nor¬ 

malized  longshore  energy  flux  defined  by 

f/  y  C  (f  d) 

E  ^  r—  AT  (16) 


which  follows  directly  from  Equation  11.  The  average  estimate  from  a  group 
of  K  constituent  spectra  within  a  class  is  defined  by 

^  i  t  '’w  0^ 

A  t.| 


When  Equation  16  is  substituted  into  Equation  17,  a  triple  summation  occurs. 
Because  each  summation  is  a  linear  combination  of  terms,  the  summation  over 
index  k  can  be  moved  inside  the  summations  over  indexes  n  and  m  with  the 
result  that  the  sum  over  k  acts  on  the  product  of  dimensionless  group  velocity 
and  constituent  spectrum.  Because  the  depth  variations  are  small  in  this  data 
set,  the  dimensionless  group  velocity  does  not  change  very  much  through  its 
dependence  on  dq)th.  Hence,  the  summation  over  k  acts  primarily  to  average 
the  constituent  spectra.  The  average  of  K  constituent  spectra  is,  by  the  defi¬ 
nition  of  Equation  7,  identical  to  the  mean  or  characteristic  spectrum  of  a 
class,  as  discussed  in  Chapter  3.  An  expression  that  should  be  very  nearly 
equivalent  to  Equation  17,  but  in  terms  of  a  characteristic  spectrum  is 

•  E  E  cos  sin  AB  Af  (18) 

11.1  ..1  ^g^ 
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where  the  overbar  on  the  subscript  indicates  the  derivation  from  a  mean  spec¬ 
trum. 


To  estimate  the  variability  of  longshore  energy  flux  estimates  from  a  set  of  K 
constituent  spectra,  a  standard  deviation  has  been  computed  for  each  class. 

For  the  K  members  of  a  class,  the  standard  deviation  is  defined  as 


p' 


E 


t«i 


K  -  1 


(19) 


Longshore  Flux  Estimates 

Figures  10  to  13  show  estimates  of  P/  from  the  unidirectional,  monochro¬ 
matic  formulation  of  Equation  14  (+’s),  the  frequency-direction  spectral  for¬ 
mulation  of  Equation  16  (x*s),  and  the  constituent  frequency-direction  spectral 
class  averages  of  Equation  17  (solid  dots)  for  all  of  the  classes 

defined  in  Chapter  2.  Each  figure  has  two  main  subplots.  Each  subplot 
contains  all  estimates  for  the  indicated  wave-height  class.  Each  subplot  has 

seven  horizontal  axes,  about  which  are  plotted  all  estimates  of  P/  as  functions 
of  9^  for  the  indicated  frequency  class  within  a  wave-height  class.  Estimates 

of  P[muji  aiid  pU*  averaged,  but  are  simply  grouped  by  class  so  that 

most  of  the  overall  data  set  is  used.  Estimates  that  fall  into  classes  for  which 
there  are  fewer  than  three  constituent  members  are  not  plotted.  Estimates  of 

pI^  are  class  averages,  and  are  plotted  along  the  direction  axis  at  dieir  class- 
center  directions. 

Because  of  the  simple  form  of  Equation  14  for  the  monochromatic,  unidi¬ 
rectional  estimates,  the  grouping  of  -H’s  that  show  these  estimates  should 
follow  a  simple  sinusoidal  pattern  with  values  of  zero  at  =  90,  0,  and 
-90  deg,  a  maximum  at  6^  =  45  deg,  and  a  minimum  at  =  -45  deg.  Ex¬ 
amination  of  Figures  10  to  13  shows  this  to  be  generally  true.  The  amplitudes 
of  the  sinusoidal  patterns  decrease  with  increasing  frequency  because  the 
frequency-dependent  group  velocity  decreases  with  increasing  frequency. 
Within  a  frequency  classification  (i.e.,  along  one  of  the  horizontal  axes  within 
a  height-specific  subplot),  there  is  some  minor  variation  in  the  values  of  the 
sinusoidal  curves  because  observed,  tidally  varying  depths  were  used  in  the 
computations,  and  becatise  each  frequency  classification  includes  a  small  range 

of  frequencies,  all  of  which  affect  the  ratio  in  Equation  14. 

In  distina  contrast,  the  sets  of  x’s  that  denote  the  frequency-direction 
spectral  estimates  form  clouds  of  points  having  considerable  variation,  and 
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Figure  11.  Longshore  energy  flux  estimates  for  wave  heights  between  0.35  and  0.84  m 
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±  0.02  Hi 


.  Longshore  energy  flux  estimates  for  wave  heights  between  1.99  and  4.70  m 


following  general  patterns  that  align  with  the  monochromatic,  unidirectional 
estimates  at  very  few  places.  The  solid  dots  that  represent  estimates  from 
characteristic  spectra  fall  within  the  clouds  of  x’s  from  the  constituent  spectra, 
as  is  expected  from  the  discussion  leading  to  Equation  18. 

These  results  mean  that  highly  resolved  spectra  give  very  different  results 
from  those  of  simple,  unidirectional,  monocfoomatic  representations  of  sea 
state.  Important  questions  that  naturally  arise  from  these  results  relate  to  the 
severity  of  the  differences  and  the  extent  to  which  corrections  can  be  made  in 
those  instances  where  well-resolved  spectra  are  not  available  and  only  charac¬ 
teristic  height,  frequency,  and  direction  are  known. 


Error  Analysis 

To  evaluate  the  differences  in  the  various  estimates  of  longshore  energy 
flux,  three  sets  of  estimators  are  compared.  First,  characteristic  spectrum 

estimates  PIj^  are  compared  to  the  means  and  standard  deviations 
from  constituent  spectra  in  their  classes  to  see  how  well  class  groupings  of 
measured  spectra  compare.  Second,  characteristic  spectrum  estimates  are 

compared  to  means  P[^  from  monochromatic,  unidirectional  estimates  in 
their  classes  to  see  if  there  is  any  general  trend  in  estimate  differences.  Final¬ 
ly,  differences  of  characteristic  spectrum  estimates  from  means  of  monochro¬ 
matic,  unidirectional  estimates  are  computed  class  by  class  to  define  empirical 
correction  functions. 

Figure  14  shows  a  correlation  of  with  class  statistical  properties  from 

constituent  frequency-direction  spectra  in  the  form  ^  ±  pU^  for  all  class¬ 
es.  Symbols  distinguish  wave  height  classes,  and,  because  high-energy  cases 
are  more  important,  solid  symbols  have  been  used  to  denote  classes  for  which 

characteristic  wave  heights  exceed  2  m.  It  is  evident  in  Figure  14  that  P'^ 

correlates  very  well  with  Pj^,  meaning  thtL.  computations  using  characteristic 
spectra  are  good  estimators  of  the  expected  or  mean  longshore  energy  fluxes 
for  all  cbsses. 

Also  evident  in  Figure  14  is  that  there  is  considerable  variation  in  energy 
flux  estimates  within  classes  as  indicated  by  the  lengths  of  the  standard  devia¬ 
tion  lines.  Standard  deviations  are  reasonably  uniform  (for  the  purposes  of 
this  study),  having  nominal  dimensionless  values  of  ±0.05  for  foe  entire 
range  of  flux  estimates.  This  result  means  that  longshore  energy  flux  esti¬ 
mates  computed  with  characteristic  spectra  will  have  an  associated,  uncorrec- 
table  uncertainty  of  this  size.  The  uncertainty  cannot  be  expressed  as  a  per¬ 
centage  because  longshore  fluxes  can  be  zero.  The  effect  of  foe  uncertainty  is 

worst  near  P/  =  0  because  varying  foe  estimate  within  foe  range  of  uncer¬ 
tainty  can  cause  flux  estimates  to  change  sign.  This  effect  is  especially  bad 
for  high-energy  cases,  several  of  which  are  located  near  foe  origin  of 
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Figure  14.  Correlation  of  P,'  estimates  from  characteristic  spectra  with 
those  of  class  averages  of  constituent  spectra 

Figure  14.  If  used  in  a  sediment  transport  model,  the  effect  of  changing  sign 
on  the  longshore  energy  flux  in  high-energy  conditions  would  result  in  a  large 
volume  of  sediment  being  moved  the  wrong  way,  which  may  be  an  unaccept¬ 
able  result  in  an  engineering  application.  This  result  suggests  that  a  more 
refined  method  of  spectral  classification  may  be  justified.  Though  the  result  is  > 
disquieting,  it  has  the  advantage  that  a  magnitude  can  be  assigned  to  the  un¬ 
certainty  in  longshore  flux  estimates  from  characteristic  spectra. 

Even  more  disquieting  are  estimates  of  longshore  energy  flux  from  Class 
averages  using  monochromatic,  unidirectional  wave  models.  Figure  IS  illus¬ 
trates  the  correlation  of  with  P^jj  for  all  classes.  Standard  deviations  for 
the  monochromatic,  unidirectional  estimates  are  not  shown  because  they  are 
very  small  relative  to  the  symbol  sizes  used  in  Figure  15.  Again,  solid  sym¬ 
bols  denote  high-energy  classes.  Though  there  is  a  vague  pattern  to  the  scat¬ 
ter  of  symbols  (monochromatic,  unidirectional  estimates  tend  to  be  larger  than 
those  from  characteristic  spectra),  there  is  not  enough  correlation  to  suggest 
that  an  accurate  correction  could  be  devised  based  on  this  representation  of  the 
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Figure  15.  Correlation  of  longshore  energy  flux  estimates  from  characteris¬ 
tic  spectra  with  those  of  class  averages  of  monochromatic, 
unidirectional  spectra 

data  set  as  a  whole.  The  differences  are  substantial.  Several  of  the  high- 
energy  symbols  fail  in  one  of  the  two  quadrants  (upper  left  and  lower  right) 
indicating  results  of  opposite  sign.  Other  high-energy  symbols  have  character¬ 
istic  spectral  values  near  zero  while  the  monochromatic,  unidirectional  esti¬ 
mates  have  values  very  different  from  zero. 

All  is  not  lost,  however.  Examination  of  Figures  10  to  13  indicates  that 
the  differences  between  the  rather  clean  pattern  of  monochromatic,  unidirec¬ 
tional  estimates  and  the  likewise  clean  pattern  of  characteristic  spectral  esti¬ 
mates  follow  a  systematic  pattern  that  is  a  function  of  the  characterizing  pa¬ 
rameter  set  .  This  pattern  is  called  a  bias,  and  is  plotted  in 

Figure  16  as  -  P',jj  as  a  function  of  direction-class  center  direction  for 
all  wave  height  and  frequency  classes,  following  the  pattern  used  in 
Figures  10  to  13.  Where  data  are  available  in  Figure  16,  it  appears  that  the 
bias  has  broad  tendencies  to  be  small  near  peak  directions  of  0  deg,  with  a 
maximum  near  45  deg,  and  a  minimum  near  -45  deg.  Subtleties  in  its  behav¬ 
ior  do  not  make  it  amenable  to  analytic  description,  however,  so  the  empirical 
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Method  of  Correction 


The  results  of  the  last  section  show  that  an  estimate  of  longshore  energy 
flux  from  a  monochromatic,  unidirectional  representation  of  a  wave  field 
defined  only  by  the  parameters  is  radier  poor  when  coiiq)ared  to 

estimates  based  on  averages  of  highly  resolved,  naturally  occurring, 
frequency-direction  spectra.  Without  correction,  engineering  computations 
based  on  monochromatic,  unidirectional  estimates  are  likely  to  be  in  serious 
error,  especially  at  times  of  high  energy.  However,  because  the  bias  function 
is  reasonably  well-behaved,  as  shown  in  Figure  16,  and  the  uncertainty  associ¬ 
ated  with  estimates  based  on  characteristic  spectra  is  quantifiable,  as  shown  in 
Figure  14,  it  is  possible  to  formulate  an  improved  estimate  that  includes  a 
measure  of  statistical  error.  The  objective  is  to  provide  a  better  estimate  of 
longshore  energy  flux  in  those  instances  when  all  that  is  known  are  three 
characteristic  wave  field  parameters. 

The  improved  estimate  consists  simply  of  a  monochromatic,  unidirectional 
estimate  from  which  bias  has  been  subtracted  and  to  which  is  affixed  an  esti¬ 
mate  of  uncertainty.  The  mathematical  expression  of  this  estimate  is,  for 
dimensionless  estimates 

cos  sin  -  bias(H^J^,e;}  ±  0.05  (20) 

igd 


where  biasiH^J^,0^  is  found  graphically  from  Figure  16,  and  indicates  a 
computational  estimate.  A  dimensional  estimate  P,  is  computed  by  multiply¬ 
ing  both  sides  of  Equation  20  with  the  normalization  factor  pga^^gd  ,  and 
noting  that  >  H^H6.  The  dimensional  result  in  terms  of  the  three  charac¬ 
teristic  parameters  and  0^,  assuming  water  density  p  and  depth  d  are 

known,  is 


A  '  ^  ^ 


cos  sin  -  bias{H^J^,9^  ±  0.05 


(21) 


Equation  21  represents  the  optimum  estimate  of  longshore  energy  flux 
based  on  three  known  parameters  because  removal  of  the  bias  converts  the 
estimate  to  the  expected  or  most  probable  value  for  each  class.  The  remaining 
uncertainty  represented  by  the  standard  deviation  estimate  can  not  be  reduced 
further  because  of  the  constraints  imposed  by  the  present  method  of  analysis. 
An  engineer  wishing  to  compute  longshore  energy  flux  can  use  Equation  21 
with  the  knowledge  that  it  gives  the  most  probable  estimate  and,  farther,  can 
determine  the  effect  of  uncertainty  in  the  estimate  on  engineering  computa¬ 
tions. 
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There  are  two  major  points  to  keep  in  mind  before  employing  Equation  21. 
First,  it  is  a  result  of  the  particular  type  of  analysis  done  in  this  r^rt.  If  the 
set  of  characterizing  parameters  is  different  (e.g.,  using  root-mean-square 
wave  height  instead  of  or  mean  wave  direction  instead  of  6^),  the  conse¬ 
quences  of  spectral  averaging,  spectral  grouping,  and  variance  estimates  are 
likely  to  be  different.  This  situation  may  require  making  additional  assunq>- 
tions  about  relationships  among  characterizing  parameters,  or  may  invalidate 
the  use  of  Equ^on  21  altogether.  Second,  the  results  shown  here  are  very 
specific  to  observations  at  the  FRF.  Application  at  sites  having  similar  forc¬ 
ing  and  boundary  conditions  will  likely  produce  meaningful  results  that,  at 
least,  will  provide  some  guidance.  Application  of  Equation  21  at  dissimilar 
sites  is  not  recommended  because  climatological  relationships  between  spectral 
parameters  and  true  spectral  shapes  is  generally  not  known.  If  such  relation¬ 
ships  are  known,  an  analysis  like  the  present  one  would  prove  useful. 

As  a  final  observation  for  this  chapter,  it  is  noted  that  the  uncertainty  in  dte 
estimate  of  longshore  energy  flux  in  Equation  21  can  be  reduced  if  additional 
parameters  are  used  to  classify  the  basic  spectra  in  the  database.  Measures  of 
such  things  as  directional  spread,  and  the  existence  and  strengths  of  multi¬ 
modal  peaks  in  the  spectra  would  result  in  smaller  variations  in  the  statistics 
of  mean  spectra,  as  discussed  in  Chapter  3.  Consequent  estimates  of  long¬ 
shore  energy  flux  would  then  vary  less  within  classes,  and  better  estimates 
could  be  computed.  The  utility  of  this  sq)proach  is  limited,  however,  because 
measures  of  these  additional  parameters  are  not  widely  done  in  field  studies, 
so  the  requisite  information  is  often  lacking.  When  such  parameters  become 
more  commonly  available,  it  would  be  most  useful  to  perform  a  more  refined 
analysis  of  the  present  data  set,  and  develop  an  estimator  for  longshore  energy 
flux  that  has  a  smaller  confidence  interval. 
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5  Radiation  Stresses 


Another  direction-sensitive  entity  used  in  modeling  and  computation  of 
nearshore  processes  is  the  radiation  stress  tensor,  which  has  three  unique 
conqmnents  that  arise  when  solutions  from  linear  wave  theory  are  ^plied  to 
the  nonlinear  terms  in  the  theoretical  equations  for  the  conservation  of  mo¬ 
mentum.  Spatial  divergences  of  radiation  stresses,  which  reach  maxima  in  the 
surf  zone,  are  primary  causes  of  wave-driven  currents.  These  divergences  are 
conq)uted  from  differences  of  radiation  stress  conq>onents  at  various  points  in 
space.  It  is  essential  that  correct  radiation  stress  values  are  used  for  proper 
computation  of  nearshore  currents,  as  well  as  processes  resulting  from  near¬ 
shore  currents,  such  as  sediment  transport. 

The  procedure  in  this  chapter  is  the  same  as  in  Ch^ter  4.  Governing 
equations  are  defined  and  estimates  of  radiation  stress  terms  by  three  methods 
are  compared.  Finally,  methods  are  suggested  for  making  optimal  estimates 
when  the  only  information  available  is  the  characteristic  parameter  set 

Definitions 

The  radiation  stress  tensor  can  be  represented  by  a  two-dimensional,  sym¬ 
metric  matrix  with  components  given  the  general  symbol  S^.  Because  of 
symmetry,  there  are  three  unique  components  5,,,  S^^  (which  equals  and 
^22 .  The  term  is  often  called  and  represents  the  onshore  flux  of  on¬ 
shore  wave  momentum.  The  term  is  often  called  and  represents  the 
onshore  flux  of  longshore  wave  momentum.  The  along^ore  flux  of  long¬ 
shore  wave  momentum  is  represented  by  or  5^.  Based  on  the  work  of 
Longuet-Higgins  and  Stewart  (1964),  Horikawa  (1988)  shows  a  compact 
formulation  for  the  radiation  stresses  for  unidirectional,  monochromatic 
waves,  and  Thornton  and  Guza  (1989)  show  the  form  of  S,,  when  the  wave 
energy  distribution  is  in  the  form  of  a  discrete  frequency-direction  spectrum 
.  Combining  these  various  expressions  results  in  a  general  equation 
for  the  radiation  stresses  in  terms  of  a  discrete  spectrum.  This  equation  is 


Chapter  5  Radiation  Stresses 


39 


(22) 


N  U 


«  E  E  ^  A/'  *./  ■  1.2 


ll»l  mm\ 


where  is  a  dimensionless  weighting  function  defined  by 

W.-A)  -  [  «/«j  ♦ «,  ]  -  i «, 


Within  the  weighting  function.  e,(dj  -  cos  6^,  ejSJ  “  6^  is  the 

Kronecker  delta  function,  which  is  equal  to  1  when  i  1°  j,  and  equal  to  0 
otherwise.  Cif^,d)  is  wave  phase  speed  from  the  dispersion  relation  of  linear 
wave  theory. 

As  with  estimates  of  longshore  energy  flux  discussed  in  Qiapter  4,  it  is 
useful  to  eliminate  variations  of  total  wave  energy  within  a  wave-height  class 
from  calculations  of  radiation  stresses,  so  that  results  will  give  a  stronger 
indication  of  directional  sensitivity.  If  both  sides  of  Equation  22  are  divided 

by  pgo^,  and  it  is  noted  that  division  by  converts  the  frequency-direction 
spectrum  to  the  normalized  form  of  Equation  6,  an  expression  for  dimension¬ 
less  radiation  stresses  is  found.  That  expression  can  be  written 

si  -  E  E  AT  IJ-IJ  <2“) 

M-l 


When  used  with  monochromatic,  unidirectional  waves.  Equation  24  reduc¬ 
es  to  a  simpler  form.  Discrete  spectra  representing  such  waves  are  zero  ev¬ 
erywhere  excqit  at/,  ■  and  ■  6^,  at  wdiich  point  S' if ^,9^  ■  \l6B£f. 

Using  these  relationships  in  Equation  24,  and  letting  5^^  rqiresent  a  dimen¬ 
sionless,  monochromatic,  unidirectional  radiation  stress  estimate,  results  in 


s' 


(25) 


To  distinguish  the  k*  such  estimate  among  a  set  of  occurrences  in  a  class, 
the  additional  subscript  k  is  attached  to  the  radiation  stress  estimate  and  to  the 
parameters  upon  which  it  depends.  The  result  is 
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C(f,j.A) 


An  average  of  K  such  estimates  within  a  class  is 


^  s' 

A  4.1 


(27) 


A  class  grouping  also  contains  K  constituent  frequency-direc- 

spectra,  of  which  the  k*  member  is  S't(f^y6^.  A  radiation  stress  estimate 
d  on  this  constituent  spectrum  is 


E  E  ^  Af 

liBl  m^l 


(28) 


Statistics  of  a  group  of  K  such  estimates  include  a  mean  Sfiji  and  standard 
deviation  of  the  set.  The  mean  is  defined  by 


■it  s', 

A.  tml 


(29) 


and  the  standard  deviation  is 


s' 


t>i 


VVB 


K  -  1 


(30) 


To  compare  the  statistics  of  Equations  29  and  30  with  the  results  of  a 
statistically  averaged  or  characteristic  spectrum,  as  defined  in  Ch^t^  3,  one 
can  compute 


SU’tt  AS  M 

n*!  M*! 


(31) 


Chapter  5  Radiation  Straasas 


41 


It  is  anticipated  that  will  correlate  well  with  because  die  nominal 
8-m  water  d^th  at  the  site  of  the  FRF  directional  wave  gauge  is  large  relative 
to  the  tide  range.  Consequently,  total  water  depth  does  not  vary  greatly,  and 
depth-dependent  functions,  like  the  ratio  of  group  velocity  to  wave  phase 
speed  in  do  not  change  much  from  observation  to  observation.  Averages 
of  radiation  str^s  estimates  from  constituent  spectra  should  thus  behave  like 
estimates  from  averaged  spectra. 

Most  of  the  discussion  of  estimate  comparisons  and  error  analysis  in  this 
ch^ter  will  be  related  to  (or  S^)  because  it  is  typical  of  radiation  stress 
components,  and  treatment  of  the  other  components  follows  the  same  pattern. 
Summary  statements  for  the  other  components  will  be  given  in  this  chapter 
based  on  supporting  figures  found  in  Appendix  C. 


Estimates  of 

For  the  computation  of  the  weighting  function  can  be  replaced  in 
Equation  24  and  following  equations  by  its  trigonometric  function  content  with 

the  delta  functions  evaluated.  This  process  results  in  the  general  form  of 
from  Equation  24  as 

*  E  E  S'tfJJ  cos  sin  AO  Af  (32) 


Note  that  this  equation  is  much  like  that  for  longshore  energy  flux,  the  dif¬ 
ference  being  that  dimensionless  spectra  are  weighted  by  ratio  of  group  veloci¬ 
ty  to  phase  speed  for  5^,  and  by  ratio  of  group  velocity  to  for  P',  . 
When  the  replacement  and  evaluation  process  is  applied  to  Equations  25  to  31, 

the  resulting  equations  relate  various  ways  to  compute  5^  and  its  statistics, 
and  have  the  same  trigonometric  and  delta  function  dependence  as  Equa¬ 
tion  32. 


S'y  Distributions 


Figures  17  to  20  show  estimates  of  Siy  from  three  methods  of  computation. 
Unidirectional,  monochromatic  estimates  from  Equation  26  are  shown  as  +’s. 
Estimates  based  on  Equation  28,  using  individual  constituem  spectra,  are 
shown  as  x’s.  Solid  dots  represent  results  from  Equation  31  using  the  charac¬ 
teristic  spectra  defined  in  Chapter  3.  Results  from  all  height,  frequency,  and 
direction  classes  are  shown  following  the  same  display  pattern  us^  in  Figures 
10  to  13. 
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igure  17.  Estimates  of  S^y'  lor  wave  heights  between  0.15  and  0.35  m 


Figure  18.  Estimates  of  for  wave  heights  between  0.35  end  0.84  m 
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Figure  19.  Estimates  of  S^y'  for  wave  heights  between  0.84  and  1.99  m 


As  with  the  results  for  longshore  energy  flux,  the  groupings  of  +  ’s  that 
represent  unidirectional,  monochromatic  estimates  follow  simple  sinusoidal 
patterns  with  rather  small  scatter.  In  terms  of  directional  dependence, 
functional  forms  have  maxima  at  45  deg,  minima  at  -45  deg,  and  zeros  at  90, 
0,  and  -90  deg.  The  amplitudes  of  the  sinusoidal  forms  vary  with  frequency 
because  of  dependence  on  the  ratio  .  Some  slight  scatter  in 

estimates  is  evident  within  classes,  and  is  due  to  modest  tidal  varia¬ 

tions  in  water  depth. 

In  contrast  with  results  for  longshore  energy  flux,  the  curves  for  mono¬ 
chromatic,  unidirectional  estimates  fall  more  nearly  in  the  centers  of  clouds  of 
x’s  that  represent  estimates  from  :onstituent  frequency-direction  spectra.  This 
result  indicates  that  there  tends  to  be  less  bias  in  estimates  based  on  mono¬ 
chromatic,  unidirectional  waves.  On  the  other  hand,  the  clouds  of  x’s  tend  to 

scatter  ov"^  broader  ranges  of  S'^ .  This  result  indicates  that  estimates  from 
individr^^  requency-direction  spectra  within  a  class  are  more  sensitive  to  the 
specif!  naoire  of  energy  distributions  within  a  spectrum  than  were  corres¬ 
ponding  estimates  of  longshore  energy  flux  in  Chapter  4. 

Solid  dots,  which  represent  S'^  estimates  from  characteristic  spectra  within 
a  class,  appear  to  correspond  well  with  centers  of  the  clouds  of  x’s,  suggest¬ 
ing  that  characteristic  spectra  can  represent  class  averages  of  results  from 
constituent  spectra  with  reasonable  verity.  The  solid  dots  are  often,  but  not 
always,  near  the  patterns  of  -f’s  of  monochromatic,  unidirectional  estimates, 

again  suggesting  a  reasonable  correlation  between  these  two  measures  of  Sl,. 


Error  analysis 

To  quantify  the  observations  made  in  the  last  part  of  this  ch^ter,  correla¬ 
tions  are  computed  for  estimates  from  characteristic  spectra  with  statistics 
from  estimates  of  constituent  spectra,  and  with  class  averages  of  unidirection¬ 
al,  monochromatic  estimates.  Also  computed  are  biases  in  unidirectional, 
monochromatic  estimates  relative  to  characteristic  spectra.  These  conq>uta- 
tions  define  a  measure  of  the  extent  of  errors  incurred  when  using  unidirec¬ 
tional,  monochromatic  estimates  of  and  provide  a  means  of  making  an 
optimum  correction. 

Figure  21  is  a  correlation  of  S'^j^  (from  Equation  31)  with  means  and 
standard  deviations  of  class  groupings  of  constituent  spectra  in  the  form 

(from  Equations  29  and  30,  respectively).  Symbols  distinguish 
wave  height  classes,  with  high-energy  >  2  m)  classes  being  shown  as 
solid  symbols.  As  anticipated.  Figure  21  shows  that  S'^jg  correlates  very  well 

with  Slyju,  meaning  that  characteristic  spectra  can  be  used  with  considerable 
success  to  estimate  the  class  average  properties  of  the  constituent  spectra. 
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Figure  21.  Correlation  of  3,^'  estimates  from  characteristic  spectra  with 
class  statistics  of  constituent  frequency-direction  spectra 

On  the  other  hand.  Figure  21  indicates  that  the  variability  of  estimates 
within  classes  of  constituent  spectra  is  significant,  as  shown  by  the  extent  of 

the  vertica!  lines  representing  There  is  no  strong  indication  that  the 

class  standard  deviations  vary  with  the  value  of  5^,  and  Figure  21  indicates  a 
nominal  value  of  +0.15  as  a  representative  standard  deviation  for  the  data  set 
as  a  whole.  This  result  shows  the  level  of  uncertainty  in  estimates  where 
characteristic  spectra  are  used.  The  uncertainty  generally  cannot  be  expressed 

as  a  percentage  because  Si,  goes  through  zero.  For  estimates  near  zero,  the 
level  of  uncertainty  is  of  considerable  importance  because  it  can  allow  the  sign 

of  to  change  in  addition  to  very  large  percentage  changes.  A  minimum 

percentage  of  uncertainty  occurs  at  extreme  values  of  5^,  and  is  about  30 
percent  for  the  data  shown  in  Figure  21. 

Figure  22  shows  the  correlation  of  5^  estimates  from  characteristic  spectra 
from  Equation  31)  with  class-averaged  unidirectional,  monochromatic 
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Figure  22.  Correlation  of  estimates  from  characteristic  spectra  with 
those  of  class  averages  of  monochromatic,  unidirectional  spectra 


estimates  from  Equation  27).  For  the  data  set  as  a  whole,  the  points 
scatter  roughly  evenly  about  the  line  of  perfect  correlation.  Close  examination 
of  Figure  22  shows  a  tendency  for  monochromatic,  unidirectional  estimates  to 

underestimate  relative  to  characteristic  spectra  at  low  characteristic  wave 

heights,  and.  more  important,  to  overestimate  5^  for  large  wave  heights. 

This  property  is  not  examined  turther  here  because  the  scatter  of  points  is 
rather  large  for  all  wave  heights,  indicating  that  there  is  not  a  generally  high 
correlation  for  the  whole  data  set.  and  that  methods  of  correction  using  the 
correlation  in  Figure  22  would  have  a  large  uncertainty. 

A  better  approach  to  making  corrections  to  monochromatic,  unidirectional 
estimates  is  to  compute  and  remove  the  biases  so  that  the  remaining  uncertain¬ 
ty  is  only  that  associated  with  variability  within  classes.  Figure  23  shows  the 

computed  biases,  which  are  defined  as  as  functions  of  direction- 

class  center  directions  for  all  wave  height  and  frequency  classes.  The  bias 
patterns  in  Figure  23  are  not  quite  as  clearly  defined  as  were  the  longshore 
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Figure  23.  Bias  in  S.y  estimates 

energy  flux  biases  of  Figure  16.  but  there  are  some  general  tendencies.  There 
is  a  tendency  for  the  biases  to  be  positive  for  low  negative  for  high 
and  often  near  zero  for  intermediate  At  low^,  there  appears  a  tenden¬ 
cy  for  a  minimum  bias  at  negative,  intermediate  6^,  and  a  maximum  bias  at 
positive,  intermediate  6^.  The  pattern  of  points  does  not  lend  itself  to  ready, 
analytical  description.  However,  Figure  23  can  be  used  as  graphical  guidance 
for  estimating  bias. 
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S,y  predictor 

To  incorporate  these  results  in  an  estimation  scheme  where  all  that  is 
known  is  the  characterizing  parameter  set  ,  one  would  compute  a 

monochromatic,  unidirectional  value,  subtraa  the  bias,  and  include  a  residual 

measure  of  uncertainty.  Using  to  represent  such  an  estimate,  the  mathe¬ 
matical  expression  is 


where  bias(H^J^,6^  is  found  from  Figure  23,  and  the  uncertainty  is  as  evi¬ 
denced  in  Figure  21.  To  create  a  dimensional  estimate  both  sides  of 
Equation  33  are  multiplied  by  the  normalizing  factor  >  pgH^l6,  so 
that  the  result  is  in  terms  of  the  three  characterizing  parameters  and 

d_.  The  estimate  takes  the  form 

P 


where  it  is  assumed 

Estimates  of 


CV.A 


that  water  density 


S«  and  S„ 


sin  «,  -  ±  0.15 

p  and  depth  d  are  known. 


(34) 


Patterns  for  the  analyses  of  the  other  two  radiation  stress  components  are 

the  same  as  for  Si,.  Results  of  these  analyses  are  simply  summarized  here, 
and,  to  save  space  in  the  body  of  this  report,  supporting  figures  can  be  found 
in  Appendix  C.  General  mathematical  forms  for  these  two  components  are 
found  by  substituting  the  appropriate  trigonometric  functions  and  evaluating 

the  delta  functions  in  Equation  23.  The  form  for  S'„  (or  S'u )  becomes 


sL 


if 


E 


S'tfJJ  M  tif 


and  the  form  for  (or  S^)  becomes 
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Equations  derived  from  Equations  25  to  31  for  monochromatic,  unidirectional 
estimates,  constituent  spectral  estimates,  characteristic  spectral  estimates,  and 
their  statistics  also  follow  form  when  the  appropriate  trigonometric  functions 
are  substituted  and  delta  functions  evaluated.  The  reader  can  readily  perform 
these  operations,  so  the  resulting  equations  are  not  listed  here. 


S„  distributions 

The  onshore  transport  of  onshore  wave  momentum  represented  by  S'„ 
tends  always  to  be  positive,  with  maxima  near  peak  directions  of  0  deg,  and 
minima  at  extreme  peak  wave  approach  angles  of  90  and  -90  deg.  This  be¬ 
havior  is  illustrated  in  Figures  (31  to  C4,  which  show  the  distributions  of 

monochromatic,  unidirectional  estimates  (+’s),  constituent  spectral 

estimates  (x’s),  and  characteristic  spectral  estimates  (solid  dots) 
for  all  cases  in  the  database.  These  figures  show  a  tendency  for  nmnochro- 
matic,  unidirectional  estimates  to  be  higher  than  the  cloud  of  points  represent¬ 
ing  spectral  estimates  near  peak  directions  of  0  deg,  and  to  be  lower  at  high 
wave  approach  angles. 

Figure  C3S  shows  the  correlation  of  characteristic  spectral  estimates 
with  the  class  means  and  standard  deviations  of  constituent  spectral  estimates 

in  the  form  i:  and  indicates  that  class  mean  values  are  well-rq)- 

resented  by  those  of  charaaeristic  spectra.  Figure  C5  also  indicates  a  general 
uncertainty  of  estimates  within  a  class  of  about  ±0.1,  whidi  can  be  expressed 

as  a  percentage  error  ranging  from  about  20  percent  for  small  to  about 

9  percent  for  large  S'„.  Figure  C6  shows  the  correlation  of  characteristic 
spectral  estimates  with  class  means  of  monochromatic,  unidirectional  esti¬ 
mates,  and  emphasizes  the  tendency  for  monochromatic,  unidirectional  esti¬ 
mates  to  be  high  for  large  S'„  and  low  for  small  S'„-  This  tendency  is  quanti¬ 
fied  class  by  class  in  Figure  C7,  which  shows  the  empirical  bias  function  of 
monochromatic,  unidirectional  class  averages  relative  to  characteristic  spectra, 

defined  as  -  S^jg. 


S„  predictor 

These  results  lead  to  a  method  of  estimation  of  5s  when  all  that  is  known 
is  the  characteristic  parameter  set  .  Designated  this  estimate  is 
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4  -  (  cos’  *  1  )  -  i  -  Wos(/C/,.»P  ±  0.1  (37) 


and  the  dimensional  form  is 

S-  ■ 


Pg 
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1(38) 


C0;.d) 


when  both  sides  of  Equation  37  are  multiplied  by  the  normalization  factor 
pgd^  -  pgHll\6. 


Syy  distributions 

The  behavior  of  S'„,  the  alongshore  transport  of  longshore  momentum,  is 

somewhat  like  the  behavior  of  the  exception  being  that  the  trigonometric 

functional  form  is  sin^  6  for  5^  instead  of  cos^  6.  This  difference  causes  5!^ 
to  have  minima  near  peak  wave  approach  angles  of  0  deg,  and  maxima  at  ap¬ 
proach  angles  of  90  and  -90  deg.  This  behavior  is  shown  in  Figures  C8  to 

Cll,  which  show  monochromatic,  unidirectional  estimates  C+’s), 

constituent  spectral  estimates  5^^  (x’s),  and  characteristic  spectral  estimates 

cases.  For  this  component  of  radiation  stress,  mono- 
c^omatic,  unidirectional  estimates  tend  to  underestimate  spectral  estimates  for 
small  peak  wave  approach  angles,  and  to  overestimate  spectral  estimates  at 
large  !q)proach  angles. 

Figure  C12  shows  the  correlation  of  characteristic  spectral  estimates 
with  the  constituent  spectral  class  means  and  standard  deviations  in  the  form 

As  for  the  other  radiation  stress  components,  this  figure 
shows  that  is  very  well  represented  by  and  that  the  class  uncertain¬ 
ty  is  reasonably  uniform  for  the  full  range  of  but  the  level  of  uncertainty 

is  smaller  than  for  S'„,  being  about  ±0.05.  Figure  C13  shows  the  correlation 
of  with  S'yy^,  and  indicates  the  tendency  for  monochromatic,  unidirec¬ 
tional  estimates  to  be  high  for  large  S'„,  and  low  for  small  5^,  much  like  the 

behavior  for  S'^  in  Figure  C6.  Figure  C14  quantifies  the  bias  -  S^jg  as 
computed  for  all  data  classes,  and  shows  the  bias  to  be  rather  well-braved. 
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Syy  predictor 


A  predictor  for  dimensionless  longshore  transport  of  longshore  wave 
momentum  evolves  from  these  results  in  the  form 


^  ±  0.05 


The  dimensional  form  becomes 

yy 
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16  ’  C(f^,d) 
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±  0.05 


when  both  sides  of  Equation  39  are  multiplied  by  the  normalization  factor 
Pg<^  -  pgfll/16. 


Discussion 

Equations  34,  38,  and  40  are  the  optimum  estimates  of  5^,  5^,  and  5^, 
respectively,  in  terms  of  the  analysis  sq)proach  used  here.  Bias  removal  con¬ 
verts  the  estimates  to  class  averages  of  constituent  frequency-direction  spectra, 
or,  very  nearly  equivalently,  estimates  from  class  characteristic  spectra.  The 
remaining  uncertainties  are  not  reducible  if  only  three  characterizing  parame¬ 
ters  are  use^  describe  the  wave  field.  It  is  likely  that  confidence  intervals 
can  be  reduc:..  if  more  parameters  are  enq)loyed  to  describe  the  constituent 
spectra.  This  last  statement  is  an  iteration  of  the  comment  made  in  Ch^ter  4, 
and  means  that  an  incomplete  description  of  the  wave  field  must  necessarily 
result  in  considerable  uncertainty  in  the  estimation  of  wave  field  properties. 

Also  iterated  here  is  the  point  that  these  results  are  specific  to  the  FRF 
and,  possibly,  similar  sites.  A  similar  site  would  be  characterized  by  adjacen¬ 
cy  to  a  broad,  shallow  continental  shelf  with  the  ocean  subject  to  wind  forcing 
by  frontal  systems  oriented  along  a  diagonal  to  the  shoreline.  While  these 
results  are  not  intended  for  general  engineering  use,  they  do  emphasize  the 
broad  truth  that  rather  significant  errors  can  occur  when  dynamically  impor¬ 
tant  properties  are  estimated  from  poorly  characterized  wave  fields. 

An  advantage  of  the  present  results  is  that  a  quantifiable  estimate  of  error 
can  be  assigned  to  radiation  stress  estimates  in  instances  where  these  results 
might  apply.  If  these  results  are  used  in  a  project  of  broader  scope,  one  can 
examine  the  consequences  of  estimation  uncertainty  on  the  project  as  a  whole. 
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and,  perhsq)s,  employ  a  compensating  safety  factor.  Because  these  results 
show  the  uncertainties  to  be  frequently  rather  large,  it  is  very  clear  that  im¬ 
proved  estimation  of  radiation  stress  conq)onents  is  possible  only  with  more 
knowledge  of  the  wave  field  than  the  three-parameter  set  that  is 

commonly  used.  Measures  of  such  things  as  directional  spread  and  multiple 
modes  would  greatly  enhance  the  spectral  classification  scheme,  and  the  result 
would  be  rather  sm^ler  error  bars  on  estimates  of  radiation  stresses. 
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6  Conclusion 


The  intent  of  this  report  has  been  to  examine  a  S-year,  6,7S9-case  database 
of  high-resolution  frequency-direction  spectra  for  three  interrelated  properties, 
all  based  on  classification  of  the  spectra  by  three  parameters:  characteristic 
wave  height,  spectral  peak  frequency,  and  spectral  peak  direction.  The  first 
property  is  the  general  distribution  of  the  observed  spectra  in  wave  height- 
frequency-direction  space.  Second  is  the  definition  of  characteristic  spectra  to 
represent  the  various  classes  of  height,  frequency,  and  direction.  Third  is  an 
examination  of  estimation  methods  for  and  variability  of  some  direction- 
sensitive  wave  properties  to  determine  the  consequences  of  using  monochro¬ 
matic,  unidirectional  waves  or  characteristic  spectra  in  lieu  of  observed  spectra 
in  the  estimates. 

The  spectra  are  classified  as  being  in  one  of  eight  geometrically  spaced 
wave  height  groups,  one  of  seven  regularly  spaced  frequency  ban^,  and  one 
of  nineteen  regularly  spaced  direction  bins.  Counting  the  number  of  spectra 
in  each  class  show^  that  low-energy  wave  fields  with  high  peak  frequencies 
tend  to  arrive  from  high  approach  angles  relative  to  shore-normal.  Those  with 
intermediate  peak  frequencies  have  ^proach  angles  closer  to  shore-normal. 
Low-energy  wave  fields  with  low  peak  frequencies  tend  to  arrive  primarily 
from  slighdy  south  of  shore-normal,  the  direction  of  nearest  {q)proach  of  Ae 
continental  ^elf  break.  High-energy  spectra  tend  to  have  low  peak  frequen¬ 
cies  and  had  peak  directions  centered  more  nearly  around  shore-normal,  as  is 
expected  due  to  the  process  of  refraction. 

None  of  these  results  are  surprising,  though  the  percentage  of  occurrence 
estimates  may  help  in  climatological  studies  of  wave  field  bdiavior.  For 
example,  the  high  percentage  of  low-energy  southerly  swell,  which  is  often 
associated  with  nearshore  sediment  accretion,  might  be  related  to  the  low 
percentage  of  high-energy  spectra  that  are  usually  credited  with  beach  erosion. 
Quantitative  models  of  subtle,  direction-dependent  behavior  in  the  dynamic 
erosion-accretion  balance  might  be  developed  when  these  results  are  used  in 
conjunction  with  bathymetric  survey  data. 

Computing  averages  of  constituent  spectra  within  classes  has  provided 
representative  or  charaaeristic  spectra  for  each  class.  Directional  energy 
distributions  of  these  characteristic  spectra  reveal  the  frequent  occurrence  of 
skewed  or  multimodal  shapes.  Such  shapes  can  have  a  strong  effect  on  near¬ 
shore  processes  driven  by  these  waves,  and  suggest  that  simple  models  of 
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continental  shelf  wave  genemion  and  transformation  may  be  inadequate  to 
account  for  these  observations. 

The  characteristic  spectra  found  in  this  study  are  expected  to  be  of  consid¬ 
erable  use  to  modelers  of  nearshore  processes,  wave  generation,  and  wave 
transformation.  For  example,  in  all  of  the  high-energy  spectral  classes,  die 
characteristic  spectra  show  distinctive  features  that  not  only  persist  dirough 
individual  storms,  but  also  ^pear  in  the  same  form  from  storm  to  storm. 
These  spectra  have  directionally  broad  low-frequency  energy  peaks  and  well- 
defined,  directionally  bimodal  high  frequency  spectral  tails.  This  characteris¬ 
tic  shape  is  different  from  what  is  expected  of  a  nearshore  spectrum,  which  is 
usually  narrow-banded  in  direction  due  to  refractive  processes.  Because  these 
spectra  occur  at  high  energy,  seeking  physical  »planations  for  them  is  an 
important  topic  for  future  research,  as  are  the  consequences  for  nearshore 
processes  being  driven  by  waves  having  this  unique  spectral  signature. 

Along  with  averages  of  constituent  spectra,  standard  deviations  of  constitu¬ 
ent  spectra  are  also  discussed.  Standard  deviation  spectra  are  found  to  be  of 
comparable  size  and  shape  to  their  corresponding  mean  spectra.  Part  of  the 
variation  is  due  to  the  rather  broad  frequency  bands  and  direction  bins  used  to 
define  spectral  classes.  This  situation  allows  constituent  spectra  of  several 
peak  frequencies  or  peak  directions  to  exist  in  die  same  class,  and  can  cause 
significant  variability  as  measured  by  the  standard  deviation  spectra. 

The  remainder  of  the  variation  is  due  to  subdeties  in  spectral  sh^e  not 
distinguished  in  a  simple,  three-parameter  classification.  Directional  side 
lobes  and  relatively  low-level  secondary  energy  modes  exist  in  some  spectra, 
but  not  in  others.  Though  secondary  relative  to  spectral  peak  parameters, 
these  features  can  have  significant  effects  on  direction-sensitive,  integral  prop¬ 
erties  of  the  spectra,  such  as  longshore  energy  fluxes  and  radiation  stresses. 
More  refined  spectral  classification  can  reduce  the  variability  due  to  these 
features  by  using  additional  diaracterizing  parameters  to  account  for  such 
things  as  directional  spread  and  secondary  modal  structure.  Though  beyond 
the  scope  of  the  present  study,  such  an  effort  is  highly  recommended  for 
future  work. 

The  final  part  of  this  study  is  a  rather  laborious,  but  quite  insightful  exami¬ 
nation  of  the  consequences  of  using  both  monochromatic,  unidirectional  waves 
and  characteristic  spectra  to  estimate  longshore  energy  flux  and  radiation 
stress  tensor  components  in  light  of  distributions  of  estimates  from  constituent 
spectra.  It  is  found  that,  without  correction,  monochromatic,  unidirectional 
waves  typically  produce  estimates  with  large  errors,  and,  in  the  cases  of  long¬ 
shore  energy  fiux  and  onshore  transport  of  longshore  momentum,  can  even 
result  in  estimates  with  the  wrong  sign.  These  are  serious  errors,  and  would 
invalidate  models  of  sediment  transport  and  nearshore  currents  that  commonly 
rely  on  these  properties. 

An  optimum  correction  is  found  for  monochromafic,  unidirectional  esti¬ 
mates  of  longshore  energy  fiux  and  radiation  stresses.  The  correction  consists 
of  adjusting  the  estimate  for  the  empirically  derived  bias  of  the  monochromat- 
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ic,  unidirectional  estimate  relative  to  either  the  mean  of  the  constituent  spec¬ 
tral  estimates  within  a  class  or,  nearly  equivalently,  the  estimate  from  the  class 
characteristic  spectrum.  The  correction  is  optimum  because  it  adjusts  the 
estimate  to  the  most  probable  value  for  the  class,  and  cannot  be  refined  further 
because  all  three  dimensions  of  the  analysis  space  (wave  height,  period,  and 
direction)  are  required  to  define  the  bias. 

Considerable  uncertainty  remains  in  the  corrected  estimates,  however, 
because  some  direction-sensitive  details  of  the  constiturat  spectra  cannot  be 
discerned  in  a  three-parameter  spectral  characterization.  The  magnitude  of  the 
remaining  uncertainty  is  very  important  because  it  still  allows  estimates  of 
longshore  energy  flux  and  onshore  transport  of  longshore  momentum  to  be  of 
opposite  sign  from  the  true  values.  This  property  of  the  results  reemphasizes 
the  need  for  characterizations  of  a  wave  field  that  are  more  refined  than  the 
conunonly  used  three  parameters.  If  wave  conditions  are  to  be  defined  with  a 
limited  set  of  parameters,  the  set  must  include  enough  information  to  distin¬ 
guish  such  things  as  directionally  broad  from  narrow  spectra,  multiple  modes, 
and  significant  directional  side  lobes  of  energy. 

Perh2q>s  the  most  important  result  of  this  study  is  the  general  conclusion 
that  characterizing  wave  conditions  at  the  FRF  with  only  three  parameters  is 
unduly  restrictive.  If  accurate  modeling  and  engineering  computations  are  to 
be  performed,  the  absolute  maximum  amount  of  available  directional  informa¬ 
tion  must  be  employed;  commonly  computed,  direction-dependent  wave  prop¬ 
erties  are  just  too  sensitive  to  directional  energy  definition  to  rely  on  low- 
resolution  wave  field  characterization.  This  conclusion  implies  that  where  a 
directional  wave  gauge  is  deployed,  low-resolution  analysis  is  likely  to  result 
in  a  poorly  described  wave  field,  and  the  problems  fou^  in  this  study  will  not 
be  solved.  Hence,  it  is  recommended  that  high-resolution  analysis  techniques 
be  used  with  directional  wave  gauges  of  any  style  in  order  to  minimize  uncer¬ 
tainties  in  direction-sensitive  wave  processes. 

Finally,  it  is  re-emphasized  that  the  specific  directional  nature  of  nearshore 
waves  discussed  in  this  rq)ort  is  true  only  at  the  FRF  site  where  measure¬ 
ments  were  made.  Some  of  the  properties  described  in  this  report  may  ^ply 
at  sites  having  forcing  and  boundary  conditions  similar  to  those  at  the  FRF, 
but  are  not  likely  to  be  generally  tq)plicable.  However,  the  broader  indica¬ 
tions  of  this  study  are  more  universal:  a  three-parameter  characterization  of 
wave  energy  spectra  is  likely  to  be  unduly  restrictive,  with  direction-sensitive 
wave  properties  suffering  from  egregiously  large  error  bars.  Considerable 
mitigation  of  this  problem  is  possible  either  with  added  parametric  descrip¬ 
tions  of  wave  spectra  or  the  direct  use  of  highly  resolved  spectra. 
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Appendix  A 

Plots  of  Characteristic  Spectra 


This  appendix  contains  illustrations  of  mean  or  characteristic  normalized 

spectra  defined  by  Equation  7  in  Chapter  3  for  all  para¬ 

metric  classes  for  which  data  are  available.  Plots  are  grouped  such  that  two 
or  more  consecutive  pages  represent  a  single  wave  height  class.  Each  page 
contains  four  direction  classes  within  its  wave  height  class,  and  all  frequency 
classes  corresponding  to  the  height  and  direction  classes.  Spaces  for  classes 
with  too  few  cases  to  average  are  included  to  preserve  the  matrix  form  of  this 
display.  To  save  space,  plot  axes  are  not  labelled.  However,  all  spectra  are 
plotted  to  the  same  scale,  and  axis  scales  are  as  shown  in  Figure  A 1 ,  which 
can  be  used  as  a  template  for  axis  notation. 


Figure  A1 .  Axis  label  template 
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Appendix  B 

Plots  of  Deviation  Spectra 


This  appendix  contains  illustrations  of  normalized  standard  deviation  spec¬ 
tra  defined  by  Equation  8  in  Chapter  3  for  all  parametric 

classes  for  which  data  are  available.  Plots  are  grouped  such  that  two  or  more 
consecutive  pages  represent  a  single  wave  height  class.  Each  page  contains 
four  direction  classes  within  its  wave  height  class,  and  all  frequency  classes 
corresponding  to  the  height  and  direction  classes.  Spaces  for  classes  with  too 
few  cases  to  average  are  included  to  preserve  the  matrix  form  of  this  display. 
To  save  space,  plot  axes  are  not  labelled.  However,  all  spectra  are  plotted  to 
the  same  scale,  and  axis  scales  are  as  shown  in  Figure  Bl,  which  can  be  used 
as  a  template  for  axis  notation. 


Figure  Bl .  Axis  label  template 
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Appendix  C 

Figures  Supporting  Chapter  5 


This  appendix  contains  figures  supporting  the  discussion  in  Ch^ter  5 
regarding  radiation  stress  tensor  conoonents  Sj^,  the  onshore  transport  of 

onshore  wave  momentum,  and  S^,  the  alongshore  transport  of  longshore 
wave  momentum.  These  figures  are  placed  here  to  reduce  the  space  required 
in  the  main  text,  and  because  the  arguments  regarding  these  radiation  stress 
components  are  rather  similar  to  arguments  given  in  Chapter  5  concerning 

radiation  stress  tensor  component  S^,  the  onshore  transport  of  longshore  wave 
momentum. 
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Figure  C3.  Estimates  of  for  wave  heights  between  0.84  and  1.99  m 
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Estimates  of  Sw'  for  wave  heights  between  1.99 


Figure  C8.  Estimates  of  Syy*  for  wave  heights  between  0.15  and  0.35  m 
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Figure  CIO.  Estimates  of  Syy'  for  wave  heights  between  0.84  end  1.99  m 
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Figure  C11.  Estimates  of  Syy*  for  wave  heights  between  1.99  and  4.70  m 


Figure  Cl  2.  Correlation  of  estimates  from  characteristic 
spectra  with  class  statistics  of  constituent 
spectra 


Figure  Cl 3.  Correlation  of  S„'  estimates  from  characteristic 
spectra  with  monochromatic,  unidirectional 
class  averages 
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Notation 


bias(H^/^,6^  A  unidirectional,  monochromatic  estimate  less  the 

corresponding  characteristic  spectral  estimate  within  a 
given  three-parameter  classification 

Wave  phase  speed  at  frequency  /,  and  water  depth  d 

Wave  group  velocity  at  frequency  and  water  depth 
d 

d  Water  depth 

d^  Water  depth  for  k*  constituent  spectrum 

Trigonometric  function 

Cj  Trigonometric  function 

e^(6)  cosine  of  an  arbitrary  angle  6 

e^(d)  sine  of  an  arbitrary  angle  B 

f  Frequency 

fd  Frequency-direction 

n*  discrete  frequency 

Characteristic  peak  frequency 

Characteristic  peak  frequency  of  k*  constituent  spec¬ 
trum 

g  Gravitational  acceleration 


Appendix  D  Notation 


H 


Wave  height  of  monochromatic,  unidirectional  wave 


Energy-based  characteristic  wave  height 
Characteristic  wave  height  of  constituent  spectrum 
i  An  index 

j  An  index 

k  An  index 

K  Upper  limit  of  index  k 

m  An  index 

M  Upper  limit  of  index  m 

mu  Monochromatic,  unidirectional 

n  An  index 

N  Upper  limit  of  index  n 

P,  Longshore  energy  flux 

P,  Corrected  estimate  of  dimensional  onshore  transport  of 
longshore  wave  momentum 

pI  Dimensionless  onshore  transport  of  longshore  wave 

momentum 

Correaed  estimate  of  dimensionless  onshore  transport 
of  longshore  wave  momentum 

P'ljg  Characteristic  spectral  estimate  of  P( 

Average  frequency-direction  spectral  estimate  of  Pj 

pU*  k*  frequency-direction  spectral  estimate  of  P/ 

Pjf,.  Standard  deviation  of  frequency-direction  spectral 

estimate  of  Pi 
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Monochromatic,  unidirectional  estimate  of  P', 


P Umt 

PUj. 

TO 

siif:) 

TO 

nej 

sUej 

s'.ifJJ 

Sy 

s' 

•’vJB 

Sgjdjt 

Syju,, 
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Average  monochromatic,  unidirectional  estimate  of  P'l 

k*  monochromatic,  unidirectional  estimate  of  pI 

Discrete  frequency  spectrum 

Mean  normalized  frequency  spectrum 

Standard  deviation  frequency  spectrum 

Discrete  direction  spectrum 

Mean  normalized  direction  spectrum 

Standard  deviation  direction  spectrum 

Discrete  frequency-direction  spectrum 

Mean  normalized  frequency-direction  spectrum 

k*  constituent  frequency-direction  spectrum 

k*  normalized  constituent  frequency-direction  spec¬ 
trum 

Standard  deviation  frequency-direction  spectrum 
General  form  of  radiation  stress  tensor  con^nent 

General  form  of  dimensionless  radiation  stress  tensor 
component 

Average  frequency-direction  spectral  estimate  of  S- 

Characteristic  spectral  estimate  of  Sh 

A:*  frequency-direction  spectral  estimate  of  S'^ 

Standard  deviation  of  frequency-direction  spectral 
estimate  of  Sj/ 
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Monochromatic,  unidirectional  estimate  of 

Average  monochromatic,  unidirectional  estimate  of  S', 

Syjmjt  monchromatic,  unidirectional  estimate  of  5^ 

Onshore  transport  of  onshore  wave  momentum 

Corrected  estimate  of  dimensional  onshore  transport  of 
onshore  wave  momentum 

Dimensionless  onshore  transport  of  onshore  wave 
momentum 

Corrected  estimate  of  dimensionless  onshore  transport 
of  onshore  wave  momentum 

^Lja  Characteristic  spectral  estimate  of 

Average  frequency-direction  spectral  estimate  of  S'„ 

k*  frequency-direction  spectral  estimate  of  Si, 

SLjd,a  Standard  deviation  of  frequency-direction  spectral 

estimate  of 


Average  monochromatic,  unidirectional  estimate  of  5^ 

sLjmj,  k*  monochromatic,  unidirectional  estimate  of  Si, 
Onshore  transport  of  longshore  wave  momentum 

Corrected  estimate  of  dimensional  onshore  transport  of 
longshore  wave  momentum 

54,  Dimensionless  onshore  transport  of  longshore  wave 
momentum 

•$4  Corrected  estimate  of  dimensionless  onshore  transport 
of  longshore  wave  momentum 

^iyjB  Charaaeristic  spectral  estimate  of 
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Average  frequency-direction  spectral  estimate  of  5^ 

frequency-direction  spectral  estimate  of  Siy 

Standard  deviation  of  frequency-direction  spectral 
estimate  of 

Average  monochromatic,  unidirectional  estimate  of 

Alongshore  transport  of  longshore  wave  momentum 

Corrected  estimate  of  dimensional  alongshore  transport 
of  longshore  wave  momentum 

Dimensionless  alongshore  transport  of  longshore  wave 
momentum 

^  Corrected  estimate  of  dimensionless  alongshore  trans¬ 

port  of  longshore  wave  momentum 

Characteristic  spectral  estimate  of 

Average  frequency-direction  spectral  estimate  of 

k*  frequency-direction  spectral  estimate  of 

Standard  deviation  of  frequency-direction  spectral 
estimate  of 

Average  monochromatic,  unidirectional  estimate  of  S'„ 
S'yy^  k*  monochromatic,  unidirectional  estimate  of 
S„  Alternate  form  for 

S'n  Alternate  form  for  5^ 

S,j  Alternate  form  for 

Sj,  Alternate  form  for  S^,  equal  to  S,j 

Alternate  form  for 
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Alternate  form  for 

Volume  of  standard  deviation  frequency-direction 
spectrum 

Dimensionless  weighting  function  in  radiation  stress 
definition 

dy  Kronecker  delta  function 

4/^  Discrete  frequency  bandwidth 

AO  Discrete  direction  arc 

p  Water  density 

a  Standard  deviation 

0^  Variance 

a]  Variance  represented  by  k*  constituent  spectrum 

0  Generic  direction  variable  representing  an  angle  mea¬ 
sured  counterclockwise  from  shore  normal 

0^  m*  discrete  direction 

6^  Charaaeristic  peak  direaion 

O^j^  Charact^tic  peak  direction  of  k*  constituent  spec¬ 

trum 
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